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The purpose of this paper is to establish a geometric scattering result for a con- 
formally invariant nonlinear wave equation on an asymptotically simple spacetime. 
The scattering operator is obtained via trace operators at null infinities. The proof 
is achieved in three steps. A priori linear estimates are obtained via an adaptation 
of the Morawetz vector field in the Schwarzschild spacetime and a method used by 
Hormander for the Goursat problem. A well-posedness result for the characteristic 
Cauchy problem on a light cone at infinity is then obtained. This requires a control 
of the nonlinearity uniform in time which comes from an estimates of the Sobolev 
constant and a decay assumption on the nonlinearity of the equation. Finally, the 
trace operators on conformal infinities are built and used to define the conformal 
scattering operator. 



Introduction 

Scattering theory is one of the most precise tools to analyze and describe the asymptotic behavior 
of solutions of evolution equations. As a consequence, it has a great importance in relativity to 
understand the influence of the geometry on the propagation of waves. Scattering in relativity 
was developed by many authors: Dimock ([17]), Dimock and Kay (dHHH) an< ^ more recently, 
Bachelot (PHEE]), Bachelot and Bachelot-Motet (0), Hafner (EZ1I2E1), Hafner-Nicolas ([3D]), 
Melnyk ([35]) and Daude HSIUFJ. 

The scattering method used by these authors relies on spectral theory: this requires the 
metric to be independent of a certain time function. It has consequently been necessary to 
develop a scattering theory which is not time dependent. As remarked by Friedlander in |24| . it 
is possible to use a conformal rescaling to study an asymptotic behavior. This method was used 
for the first time by Baez-Segal-Zhou in [5] for the wave equation on the Minkowski spacetime. 
Their method consisted in embedding conformally the Minkowski spacetime in a bigger compact 
manifold. 

The conformal compactification of a spacetime was first introduced by Penrose in the sixties in 
|40) to describe the asymptotic behavior of solutions of the Dirac equations. A boundary, which 
represents in some way the infinity for causal curves, is added to the manifold. This boundary 
is divided into two connex components and J?~ . When the spacetime satisfies Einstein 
equations (with no cosmological constant), <y + and Jf~ are light cones from two singularities 
i + and i~ . The asymptotic behavior can then be obtained by considering the traces on these 
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hypersurfaces of a solution of the conformal wave equation (more precisely the scattering operator 
is obtained from the trace operators on and J ! ~). Asymptotically simple curved spacetimes, 
that is to say spacetimes admitting a conformal compactification, with specifiable regularity at 
i + and were constructed by Chrusciel-Delay 112) ). Corvino ([13]) and Corvino-Schoen 
(|14j). Mason and Nicolas successfully adapted the method of Baez-Segal-Zhou in the linear 
case for the Dirac and Maxwell equation in [3J] on this curved background. They also obtained 
a complete peeling result for the wave equation on a Schwarzschild background in [35]. 

This paper presents the construction of a conformal scattering operator for the conformally 
invariant defocusing cubic wave equation: 

V Q V a $ + 6$ 3 = 

on the asymptotically simple spacetime obtained by Chrusciel-Delay and Corvino-Schoen. Our 
construction relies on vector fields methods which were previously used to obtain the well- 
posedness of the Cauchy problem equation (see the result of Cagnac-Choquet-Bruhat in [5]) to 
obtain energy estimates. The choice which is made here is for the vector field is the same as the 
one made in [35J and the techniques are essentially the same as in |34| . 

A specific method is used to handle the singularity in i + : the characteristic hypersurface is 
described as the graph of a function. This method was introduced by Hormander in [32] and 
generalized in [3SJ [35] by Nicolas to establish the well-posedness of the characteristic Cauchy 
problem. 

The main obstacle and difference with the linear case are the necessity to obtain uniform 
estimates of the non linearity. This requires to obtain Sobolev embeddings from H 1 into L 6 over 
a spacelike foliation transverse to a light cone up to the vertex of this cone. This is achieved 
by considering results concerning the norm of this embeddings given by Miiller zum Hagen and 
Dossa (see [Till [37] an< A EH]) : studying the characteristic Cauchy problem for nonlinear equations, 
they noticed that the norm of this embedding is not bounded when the foliation reaches the 
vertex of the light cone and gave a precise estimate of the speed of the blow up. Their result is 
adapted here to the conformal infinity. The blow up of the norm of the Sobolev embedding is 
compensated by a decay assumption on the function b. 

The paper is organized as follows: 

• the first section introduces the geometrical and analytical background: the spacetimes 
obtained by Corvino-Schoen and Chrusciel-Delay are precisely defined and the function 
space on the characteristic hypersurface at infinity is given. 

• The a priori estimates are derived in section 2: these estimates are established in three 
specific subsets of M: a neighborhood of i° where the estimates come from the asymptotic 
behavior of the chosen vector field (the Morawetz vector field), a neighborhood of i + where 
the estimates are established by following the method developed by Hormander, and finally 
in a neighborhood of a Cauchy hypersurface. The techniques consist essentially in the use 
of Gronwall lemma and Stokes theorem. 

• Section 3 is devoted to the well-posedness of the Cauchy problem. The proof is made 
as follows: estimates on the propagator of the cubic wave equation are established from 
a control of the norm of Sobolev embeddings, in the spirit of the work of Miiller zum 
Hagen and Dossa. The second step consists in establishing a local existence theorem for 
the characteristic Cauchy problem up to a spacelike hypersurface close to the conformal 
infinity. This is done by slowing down the propagation of waves, reducing the characteristic 
Cauchy problem to a Cauchy problem. A complete solution, up to a reference Cauchy 
surface, is then obtained by gluing the local solution and a solution of the Cauchy problem 
on the well chosen hypersurface close to conformal infinity. 
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• The appendix contains, on the one hand, a certain number of calculations useful for the a 
priori estimates and, on the other hand, a proof of the characteristic Cauchy problem for 
small data; this proof is considered as interesting because of its simplicity. 
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Let us consider on M the following nonlinear wave equation: 

V a V a $ + 6$ 3 = 

We assume that: 

a. b is positive; 

b. b admits a continuous extension to M such that b vanishes at J?\ 

c. b satisfies: there exists a constant c such that, uniformly on M: 

3c, \f a V a b\ < cb. 

d. b satisfies the following decay assumption: there exist a uniformly spacelike foliation 
(S r ) rg [o.fl] of a neighborhood i + and a constant c, which does not depend on r such 
that 

Vre(0ijR]i M^) < c . 

The same decay assumption holds in i~ . 
These assumptions will be referred to as assumptions |SJ [bj and [d] 
Remark 0.1. a. The positivity of b corresponds to the defocusing case. 

b. Since M is compact, the differential inequality satisfied by b does in fact not impose another 
specific asymptotic behaviour than the fact that T a \7 a (j) decrease and vanishes at the same 
rate of b. This assumption, essentially technical, is required to insure the global well- 
posedness of the Cauchy problem accordingly to the Cagnac-Choquet-Bruhat theorem (see 
theorem \TJl\) . 

c. The purpose of the decay assumption on b is to compensate the blow-up of the Sobolev 
constant at the vertex of the light cones J? ± . 

The following notations will be used: 

• we will note: 

<t> < ip 

where </> and tp are two functions over U, a subset of M, whenever there exists a constant C, 
depending only on the geometry, the vector T a , the Killing form V^T^ and the function 
b, such as: 

ip < C<p on U. 

If ip and <f> both satisfy: 

<j) < i\) and ip < (j), 
we say that 4> an d ip are equivalent and note: 

• The geometric notations are the following: 

— The quantities with * are geometric quantities related to the unphysical metric. 

— fi[g] is the volume form associated with the metric g. 

— If v is a form over M, then *v is its Hodge dual. If v is a 1-form and V the vector 
field associated with v via the metric g, then: 

*v = Vjfi[g] or -kV a = V a jfj,[g] 

where jfi[g] is the contraction with the volume form fj,[g] 

— i£ is the restriction to the submanifold S. The pull-forward of a form v on M over 
the tangent space to E is denoted by i^,{v)- 
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1 Functional and geometric preliminaries 



We present in this section the geometric and analytic background. A specific care is brought to 
the structure at null infinity and the definition of function spaces on this structure. 

1.1 Geometric framework 

The geometric framework is based on the results of Corvino-Schoen ( |141 [T5] ) and Chrusciel- 
Delay([TTl[T2]). They gave a construction of asymptotically simple spacetimes satisfying Einstein 
equations with specifiable regularity at null and timelike infinities. 

1.1.1 Asymptotic simplicity 

The notion of asymptotically simple spacetimes was introduced by Penrose as a general model for 
asymptotically flat Einstein spacetimes and their conformal compactification (see |41| definition 
9.6.11): 

Definition 1.1. A smooth Lorentzian manifold M satisfying Einstein equations is said to be 
(C ) asymptotically simple if there exists a smooth Lorentzian manifold M with boundary, a 
metric g and a conformal factor il such that: 

a. M is the interior of M ; 

b. g = Q 2 g in M; 

c. O and g and are C k on M; 

d. O is positive in M ; 51 vanishes at the boundary J? of M and <Kl does not vanish at J 1 ' ; 

e. every null geodesic in M acquires a past and future end-point in 

We assume that the boundary J? is C 2 (which is sufficient for this work). It is known 
that this boundary is a null hypersurface (that it is to say that the restriction of the metric 
to J 1 is degenerate) provided that the cosmological constant is zero. Furthermore, ^ has two 
connected components and consisting of, respectively, the future and past endpoints of 
null geodesies. ^ + and are both diffeomorphic to K x S 2 . 

The manifold (M, g) is usually referred to as the physical spacetime and its compactification 
is referred to as the unphysical spacetime. In order to remain consistent with this notation all 
along this paper, the quantities associated with the unphysical metric are denoted with a 

1.1.2 Global hyperbolicity 

An important assumption in the context of the Cauchy problem for a wave equation is the 
possibility to write the equation as an evolution partial differential equation. This is usually 
achieved by requiring that the manifold M is globally hyperbolic: 

Definition 1.2. A Lorentzian manifold (M,g) is said to be globally hyperbolic if, and only if, 
there exists in M a global Cauchy hypersurface, i.e. a spacelike hypersurface such that any 
inextendible timelike curve intersects this surface exactly once. 

A useful consequence of this is the existence of a time function on M and the parallelizability 
of M, that is to say the existence of a global section of the principal bundle of orthonormal frames 
(see the work of Geroch in [25l[26] and Bernal-Sanchez in [7]). 

In the case of an asymptotically simple manifold (M,g), this property extends of course to 
the manifold (M,g). 
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1.1.3 Corvino-Schoen/Chrusciel-Delay spacetimes 

We can then introduce the spacetimes obtained by Corvino-Schoen and Chrusciel-Delay: 
Definition 1.3. A spacetime M is of Chrusciel-Delay /Corvino-Schoen type if: 

a. AI is asymptotically simple; 

b. M is globally hyperbolic; let Eo be a spacelike Cauchy hypersurface; M is then diffeomorphic 
toRx E ; 

c. M can be completed into a compact manifold by adding three points i , i + and i~ such 
that ^ + and are respectively the past null and future null cones from i + and i~~ and 
i is the conformal infinity of the spacelike hypersurface Eo for the metric g; 

d. there exists a compact set K in Eo such that (RxEo\-ftT, g) is isometric to (Kx]ro, +oo[x§ 2 , gs), 
where g$ is the Schwarzschild metric with mass m and r$ > 2m; 

e. there exists a neighborhood of i + such that the metric g is obtained in this neighborhood 
as the restriction of a smooth (C 2 ) Lorentzian metric of an extension of M in the given 
neighborhood of i + . The same property holds in i~ . 

Remark 1.4. a. The result of Corvino-Schoen/Chrusciel-Delay states that the metric is iso- 
metric to the Kerr metric outside a compact set; we restrict ourself to a Schwarzschild 
metric. 

b. The extension of the manifold M in the neighborhood of i + was used by Mason-Nicolas 
(see \34\ \35^ ). The point i + remains singular in M but it is nonetheless possible, because of 
the existence of this extension, to consider geometric data in i + (metric, exponential map, 
connection, curvature) as being the one obtained from the Lorentzian manifold extending 
(M , g) in a neighborhood of this point. 

c. The point i is a "real" singularity of M; nonetheless, the geometry is completely known in 
its neighborhood. This singularity is the main problem we have to deal with in this paper, 
being an obstacle to global estimates and Sobolev embeddings for instance. 

The last part of this section is devoted to the geometric description around the point i°, that 
is to say in the Schwarzschildean part of the manifold. 

We consider here a neighborhood O in M of i where the metric g is the Schwarzschild 
metric. This metric, in the Schwarzschild coordinates (i,r, lo^i), can be written: 

g = F(r)dt 2 - T^rdr 2 - r 2 d 2 oj S 2 
F(r) 

where: 

F(r) = 1 

r 

with m a positive constant. Introducing the new coordinates: 

r* = r + 2m log(r — 2m), u = t — r* and R = — 

r 



the metric can then be expressed as: 

2 1 

— -dudR ■„ 

R 2 R 2 



g = {l- 2mR)d 2 u - —dudR - ^d 2 ^- 



where d 2 oj§2 stands for the standard volume form on the 2-sphere, which can be written in polar 
coordinates {9, (f): 

d 2 w«2 = sin6»d6» A dd>. 



Its inverse is: 

9' 1 = J^ d ^ 2mR)d 2 R - -^99'. 

The part of the Cauchy hypersurface So is given by the equation {t = 0} in these coordinates. 
In the neighborhood O, the conformal factor is chosen to be: 

Q = R 

and is extended smoothly in M\0. The conformal metric is then: 

g = R 2 (l - 2mR)d 2 u - 2dudR - d 2 w S 2. 

and its inverse: 

g- 1 = -2d u d R - R 2 {1 - 2mR)d 2 R - d& . 

The point i is given in these coordinates by u = — oo, R = 0. 

The description of the geometry around i is completed by the following lemma (lemma A.l 
in [35]): 

Lemma 1.5. Let e > 0. There exists uq < 0, |uo| large enough, such that the following decay 
estimates in the coordinate (it, r, 9, ip) hold: 

r < r* < (1 + e)r, 1 < Rr* < 1 + e, < R\u\ < 1 + e, 

1 - e < 1 - 2mR < 1, < s = M < 1 

Proof. The proof of this lemma is straightforward: it only consists in writing simultaneously 
the asymptotic behavior or the continuity over M of each of the coordinates involved in the 
lemma. □ 

Remark 1.6. a. As mentioned in the introduction, we choose to work in the neighborhood 
of i° with the Morawetz vector T a defined by: 

f a = u 2 d u - 2(1 + uR)d R . 

The squared norm of this vector is: 

g ab f a f b = it 2 (4(l + uR) + u 2 R 2 (l - 2mR)). 

This polynomial in uR vanishes at: 



1 ± y/2mR -2 
2- 



1 - 2mR l T V2rni? 

so that, if R is chosen small enough, these roots are arbitrarily close to —2. Let e be a 
positive number chosen such that the inequalities in lemma \775\ hold for a given Uq. The 
largest zero of this polynomial satisfies: 

-2 -2 
l + V2rn~R ~ 1 + y/e' 

Choosing e such that 

-2 

< -1-e, 



1 + 

the norm of T a is then uniformly controlled by: 

g ab f a f b = u 2 (4(l + uR) + u 2 R 2 (l - 2mR)) > 4u 2 Q e 
on the domain 0+ = {(u, R,u>§?)\u < uq} PI J + (So); ^ + (^o) being the future o/So- 
b. Another criterion, given in proposition\2.4\ will be required to define e. 
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1.2 Analytical requirements 

We introduce in this section the technical and analytical tools which are required to the descrip- 
tion of the solution for the wave equation. 

1.2.1 Conformal wave equation 

We recall here how the problem on the physical space time is transformed into a problem on the 
unphysical spacetime. This is based on the classic transformation of the wave d'Alembertian 
operator: 

Lemma 1.7. Let M be a Lorentzian manifold endowed with the conformal metrics g and g = Q 2 g 
where Q is a conformal factor in C 2 (M). The connections associated with g and g are denoted 
by V and V respectively. 

Then, for any smooth function <f> on M , the following equality holds: 

v„v a </> + ^Scaigcj) = ir 3 ^v a v a (n~V) + ^Scai g (n-\ 

where Scal g and Scal g are the scalar curvatures associated with g and g respectively. 

Assuming that we are working on a vacuum spacetime, for which the scalar curvature van- 
ishes, the equation becomes: 

v a v> = rr 3 (V a v Q (jrV) + ^s c &\ g (frV)J (1.1) 

We obtain in particular the useful formula: 

n 3 W a W a n = -Scalg. (1.2) 
Let us now consider the Cauchy problem on the physical spacetime M: 

U4> + b(/) 3 = 

0| Eo ^eC °°(Sc) . (1.3) 
^Mso = £eC °°(£ ) 

Using this transformation, this Cauchy problem is transformed into a Cauchy problem on the 
unphysical spacetime M as follows: 

Lemma 1.8. The function <j) is a solution of problem (|1.3p if and only if, the function 

ip = rrV 

is solution of the problem on So : 

' hip + \Scal g ^ + bi' 3 = 

, fa, M So = n - Cr°v«n)£) G C °°(S ) 

Remark 1.9. a. Because of the finite propagation speed, since the data on the physical space- 
time are smooth with compact support on So, the data on the unphysical spacetime are 
smooth with compact support in So- 

b. Another consequence of the finite propagation speed is that, since the data remain with 
compact support in So, we do not have to deal with the singularity in i° . 
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c. Conversely, it is possible to start with a Cauchy problem on Eo in the unphysical spacetime 
and obtain a Cauchy problem on the physical spacetime: starting with the Cauchy problem 
on M: 

Dip + lScal s ip + bip s = 
^| So =0eC o °°(£ o ) 
T a V a v4 o = |eC °°(So) 

then (j) = fltjj satisfies the Cauchy problem: 

□0 + bil.' 3 = 
< ^| Eo =f^GCg°(E ) 

, ^ av ^| So = n£+(f a \7 a n)d e c °°(s ) 

1.2.2 Function spaces 

The purpose of this section is to give a precise description of the Sobolev spaces which are used 
in the present paper. Two problems are encountered in this section: the first one consists in 
adapting the definition of Sobolev spaces on a null hypersurface and the second is the difficulty 
coming from the singularity at i°. This difficulty has two aspects: the necessity to adapt the 
definition of the Sobolev space to the singularity: this is solved using weighted Sobolev spaces 
on J? . The other aspect of this singularity is encountered in section l3.1.2l when trying to obtain 
uniform Sobolev estimates of the non linearity. 

We recall the definition of a Sobolev space on a uniformly spacelike hypersurface E: for a 
smooth function u on E, consider the norm, when the integral exists: 




D k u\\ 2 h dfi[h 



where h is the restriction of the metric g and D is the restriction of the connection V to E. 
Definition 1.10. The completion of the space: 

{ue C°°(£)| ||u||p < +00} 

in the norm || * \ \ p is denoted by iJ p (E). 

When E is a compact spacelike hypersurface with boundary, the completion of the space of smooth 
functions with compact support in the interior of 'E in the norm || * \ \ p is denoted by Hq(E). 

Remark 1.11. It is known that, when working on a Riemannian closed manifold, the Sobolev 
spaces are independent of the choice of the metric. Nonetheless, this fact is not true any more 
when working with a weakly spacelike hypersurface, as we are about to see. Arbitrary choices are 
made for their definitions. 

Because of the degeneracy of the metric, it is not possible to define on geometric quan- 
tities that only depend on the metric g. Two solutions can be provided: 

• lifting the metric from Eo to J^ + ; 

• adding geometric information on by using the uniformly timelike vector field T a . 

Following [33] and using Geroch-Held-Penrose formalism, J? + is endowed with a basis 
(l a ,n a ,e1,el) such that: 
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• ^ a and n° are two future directed null vectors; n a is tangent to J^ + ; they satisfy: 

/ + n = 1 ; 
in the neighborhood of i , they are chosen to be: 

l a = -2d R and n a = u 2 d u 

• the set {l a ,n a } is completed by two vectors (ej^ef) orthogonal to {Z a ,n a }, orthogonal to 
each other and normalized. 

,y + is then endowed with the volume form i\ + (l a j(i[g]). The following norm is defined on , 
for u a smooth function with compact support which does not contain i or i + : 

lHHS + )= l \ : g cdfcfd 

where V§2U stands for the derivatives with respect to (e^,ef). 

Following chapter 5.4.3 of Friedlander's book (|23J), the Sobolev space H 1 on ,JP + is finally 
defined: 

Definition 1.12. Let M be an extension of M behind i + and consider the function space T>(^ + ) 
on .y + obtained as the trace of smooth functions with compact support in M which does not 
contain i . The weighted Sobolev space H 1 (J' + ) is defined as the completion of the space T)(^ + ) 
in the norm || * ||fr 1 (^'+)- 

Since the volume form is written on the Schwarzschildean part of M as, using polar coordi- 
nates: 

fi\g] = sm(9)du A dr A dO A dip, 
this norm is written in this region of the manifold as: 

= 2 J (ju 2 {d u <P) 2 + |V S 20| 2 + <f> 2 ^ dudw S 2 

since, on J? + , 

n a = u 2 d u and g cd T c T d = 4u 2 

The metric at i + is obtained as the restriction of a smooth metric of an extension of M 
beyond . As a consequence, a trace theorem could give another way, more intrinsic, to define 
H 1 (J? + n O) where O is a bounded open set around i + . Another way to obtain the fact that 
the point i + does not matter in the definition of the Sobolev space over J? + is to notice the 
following property: 

Proposition 1.13. The set of smooth functions with compact support on , whose support 
does not contain i + , is dense in H 1 (J? + ). 

Proof. The method of the proof relies on the construction of an identification between Hq(Eq) 
and H 1 {,# + ). This identification is brought by Hormander in [32] and is obtained as follows. 

Let f be a smooth time function in the future of Eo in M . This time function gives rise to 
local coordinates, where t is the time coordinate. We denote by dt the vector field associated 
with this coordinate. The flow associated with this vector field is denoted ty t . 

For x in Eo, let 4>{x) be the time at which the curve ^t{x) hits J? + and consider the 
application defined by: 

£ 6 C °°(E ) — ► (y G y + h+ £ (*_*(„)(!,))) 
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This application has value in C^°(J? + ) since the future of a compact set in Eo has compact 
support in . Furthermore, this application can easily be inverted: 

£ e C °°(^ + ) h_> (x e .y + ^ € (#_ 0(a) (aO)) e C °°(E ) 

and can consequently be used to define on J? + a Sobolev space by pushing forward the iJ 1 -norm 
on So. The Sobolev spaces which are obtained are then equivalent on H 1 (J? + ) since the norms 
are equivalent on any compact set of . 

To prove the proposition, it is then sufficient to prove that the smooth functions with compact 
support in Eo which does not contain the preimage of i + by the flow associated with the time 
function t. Since Eo has no topology, it is sufficient to establish the following lemma: 

Lemma 1.14. Let us consider the set of smooth functions defined in B(0, 1) C M 3 with support 
which does not contain 0. Then this set is dense in i? 1 (B(0, 1)). 

Proof. The proof is given in appendix IA.ll □ 
Using this lemma in the neighborhood of the preimage of i + immediatly gives the result. □ 

1.3 Cauchy problem 

A well-posedness result for the Cauchy problem in our framework is now stated. It is based on 
a result of Cagnac and Choquet-Bruhats in [B] (and see also [S], appendix III for a survey on 
the wave equation, and appendix III chapter 5 for our problem). 

The geometric framework for this well-posedness theorem is the following (definition 11.8 in 

my- 

Definition 1.15. A spacetime (M,g) is said to be regularly sliced if there exists a smooth 3- 
manifold E endowed with coordinates (x 1 ) and an interval I o/R such that M is diffeomorphic 
to I x E and the metric g can be written: 

g = N 2 dt 2 - g v] {dx l + f3 l dt). 

and its coefficients satisfy: 

a. the lapse function N is bounded above and below by two positives constants: 

3(c,C),0 < c < N < C; 

b. for t in I, the 3-dimensional Riemannian manifolds ^{t} X M, gt,i 
plete and the metrics g t are bounded below by a metric h i.e.: 

:\ I y.. h , , I ' I ' • .<//.,,. I T'': 

c. and, finally, the norm for the metric gt of the vector f3 is uniformly bounded on M . 

Remark 1.16. a. This hypothesis implies that the spacetime is globally hyperbolic (theorem 
11.10 in JBj). 

b. The asymptotically simple spacetime and its compacification which we are working with do 
not satisfy this property. 

The following theorem, obtained by Cagnac and Choquet-Bruhat in [5], gives existence and 
uniqueness of solutions to the Cauchy problem for a cubic wave equation: 



9 1] 

{t}xM 
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Theorem 1.17 (Cauchy problem for a nonlinear wave equation). Let us consider the Cauchy 
problem on the regularly sliced manifold (AI = lx E, g): 

Dcf>+ ±Scal a 4> + b(t> 3 = 
where the function b satisfies the assumptions^ and\Q 

Then this problem admits a unique global solution on AI in C°(R, H l (Yi)) PI C 1 (IR, L 2 (S)). 

Proof. As already noted (see remark |1.16|) , this theorem cannot of course be applied directly to 
the compactification of M because of the geometry in the Schwarzschildean part. This problem 
can be solved using an extension of AI constructed as follows: 

a. Let 9 and 9 be two functions respectively in i? 1 (Eo) and L 2 (T,q) with compact support in 
the interior of So. Let if be a compact subset of Eo containing the supports of 9 and 9. 

b. Let t be a time function on M such that Eo is given by {t = 0}; the associated foliation 
is denoted by (E t -) for t G [0,T]; we assume that the gradient of this time function is 
uniformly timelike for the metric g. 

c. The manifold (J + (K), g) is a 4-dimensional Lorentzian manifold with boundary. This 
boundary is constituted of K, the light cone from K, C + (K), and the part of J? U {i + } 
in the future of K. Since (AI,g) is extendible smoothly in the neighborhood of i + , there 
exists a smooth extension of (J + (K),g) into a 4-dimensional Lorentzian manifold (M,g), 
depending on the support of K such that the manifold M is diffeomorphic to [0,T] x E, 
where E is topologically equivalent to § 3 and such that the foliation (E 4 ~) is extended into 
the uniformly spacelike foliation ({i} x E). 

Since AI is compact, conditions 1, 2 and 3 of definition 1 1 . 1 5l are satisfied. As a consequence, 
using theorem 11.171 we obtain the well-posedness in C°(IR, if 1 (Eo)) of the Cauchy problem for 

n<f>+ iScal^ + frc/) 3 = 

^l{o}xs = ^ ^ ^H^) w ith compact support in K 
<9 t 0|| O j xS = 9 G i 2 (E) with compact support in K 

As a consequence, we obtain by restriction to J + (K) and on M well posedness of the Cauchy 
problem for data in Hq{T, ) x L 2 (E ). □ 

Remark 1.18. The extension of the solution of the wave equation to a cylinder was also used 
by Mason-Nicolas in \"33§ (proposition 6.1) to obtain energy estimates. 

2 A priori estimates 

The purpose of this section is to establish a priori estimates for solutions of the wave equation, in 
the sense that it is possible to control the energy on ,f + by the energy on Eo and reciprocally. 
These a priori estimates will be used in the next section to establish the continuity of the 
conformal wave operator and obtain its domain of definition and the existence of trace operators. 

Let us consider in this section a smooth solution with compactly supported data of the 
problem: 

Q/>+ iscal^ + 60 3 = (2.1) 
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and the associated stress-energy tensor: 

Tab = V a <t>S7 b <f> + g ab f --V c 0V c </> + — + &— 

The contraction of this tensor with T a , T a T ab , is called "energy 3-form"; it satisfies an "approx- 
imate conservation law": 

Lemma 2.1. The derivative of the energy 3-form satisfies: 

V Q (f»T ab ) = (V Q T ")) T ab ^Scal^j cf>f Q V a + T a V a b^. 

This derivative will be designated as "the error term" since it arises in the volume term when 
applying Stokes theorem. 

A quantity which is equivalent to the integral: 

i*{*T a T ab ) 

on the hypersurface S is called an energy and is denoted by E(S). Global energies on Eo and 
are piecewise defined in propositions 12.41 12. 151 and 12. 201 The purpose of this multiple definitions 
is to simplify the comparison with other quantities. 

2.1 Estimates in the neighborhood of i° 

The purpose of this section is to obtain a priori estimates for the energy associated with the 
energy 3-form on and Eo- using the fact that the geometry is known almost completely. 

The estimates which are obtained in this section are close to the ones obtained for the linear 
wave equation by Mason-Nicolas in [35] ■ These estimates are based on two main tools: 

• an explicit control of the decay of the physical metric in the neighborhood of i° 

• and the use of Gronwall lemma. 

We define, in 0+ o = {t > 0,u < uq}, the following hypersurfaces, for uo given in R: 

• S Uo = {u = uq}, a null hypersurface transverse to 



u ■ 



• Eq° = Eo D {uq > u}, the part of the initial data surface Eo beyond S., 

• j?+ g = 0+ n J+, the part of J + beyond 5„ ; 

• "K s = S7+ Q n {u = — sr*}, for s in [0,1], a foliation of f2+ by spacelike hypersurfaces 
accumulating on J^. 

The volume form associated with g in the coordinates (i?, it,w§2) is then: 

fj,\g] = du A dR A d 2 w S 2 . (2.2) 

Finally, we consider the approximate conformal Killing vector field T a : 

f a = u 2 d u - 2(1 + uR)d R . 

Remark 2.2. a. This vector field is timelike for the unphysical metric and, as such, is trans- 
verse to J?. More precisely, it is uniformly timelike in a neighborhood of i (see remark 
\1.6] for the choice of the neighborhood). 
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Figure 1 : Neighborhood of i° 

b. Its expression is derived from the so-called Morawetz vector field in the Minkowski space 
and was previously used to obtain pointwise estimates in the flat case. 

The strategy of the proof in this section is the following: 

a. writing an explicit description of the hypersurfaces S Uo and ; 

b. proving energy equalities in both ways for J -kT a T a t„ using Stokes theorem between the 
hypersurfaces and 'K" 



("HO . 



c. determining an energy E{"K^ a ) equivalent to J -kT a T a ] J from the decay of the metric g; 

d. obtaining an integral inequality for E(3i^°) to apply the Gronwall lemma; 

e. starting from point |bl doing the same work from Stokes theorem applied between .J?^[ and 



2.1.1 Geometric description 

This section is devoted to the description of the energy associated with the nonlinear wave 
equation in the neighborhood of i°. 

Proposition 2.3. The energy 3-form, written in the coordinates (R,u,9,if>), is given by: 



*T a T ab = 



u 2 (d u (/)) 2 + R 2 (l - 2mR) (u 2 d R (j)d u (j) - (1 + uR)(d R <£) 2 ) 



4 



sin(0) du A 60 A dip 



X - ((2 + uR) 2 - 2mR 3 u 2 ) (d B ^ 2 ) + u 2 Qp 



— + b— 

2 4 



sin(9)dRA dd A dijj 



+ sin(6») [u 2 d u <j) - 2(1 + uR)d R <j>] (-dg(j)du A dR A # + d^4>du A dR A 
The restriction of the energy 3-form can be written: 
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• to Vi , 



i*K a {*T a T ab ) = {u 2 {d u c^f + i? 2 (l - 2mR)u 2 d B ,<j)d u (f ) 

?, x /(2 + ui?) 2 mu 2 R 3 . ,\ .„ 

+i? 2 1 - 2mi? i * (1 + uR) (d R <j) 2 ) 

\ 2s s J 

-2mi2) , 2(1 + ui?) \ (l|v^| + ^ + &^^sin(0)duAd0A#; 



s v ' \2' a r ' 2 4 



i* Su (*f a T ab ) = I ^ ((2 + ui?) 2 - 2mi?V) (c^ 2 ) 
+m 2 Q|V S 20| 2 + ^- +b^) J sin(0)dflA d9 A #; 



i^+{-kf a T ab ) = {u 2 {d u (j)) 2 + |V S 2(/)| + 2 + 6y^) sm(8)duA d0 A dtp. 

Proof. The proof of this proposition is a straightforward calculation and is given in the appendix 
(see proposition IA.2|) . □ 

Proposition 2.4. There exists uq, such that the following energy estimates holds on ^K s in f2+ : 

A2 j,4^ 



I 



^ (*T a T Qb ) » (V(d u 0) 2 + ^(drf + |V S2 0| 2 + t. + b^jduA dLo S 2 



Proof. The proof of this proposition can be found in the appendix (see proposition IA.4jl . □ 

uo is now fixed, being equal to the uq associated with the e which ensures that the energy 
equivalence established in proposition 12.41 holds. The neighborhood of i° where the energy 
estimates are relevant is then f2+ . 

2.1.2 Energy estimates near i° 

The energy estimates are established between Sq 0> j an< ^ ^u > writing a Stokes theorem 
between ^K s , if?* = {(u, R, u>s 2 )\ u = u = uq,u < —sr*} and . 
The first step consists in evaluating the error term: 

Lemma 2.5. The error is given by: 
V a (f b T ab ) = AmR 2 (3 + uR) {d R (j)) 2 + (1 - YlmR) 4> (u 2 d u <j) - 2(1 + uR)d R <j)) + f a V a b^-. 



Proof. The proof of this lemma is essentially based on a computation (See the complete calcu- 
lation in the appendix in proposition IA.5[) □ 

Remark 2.6. As noticed in ]35 } j . one obstacle to the use of the parameter s for the foliation is 
the fact that this parametrization in s is not smooth in the sense that (r*)^ 1 is not a smooth 
function of R at R = 0. 
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In order to avoid this singularity, the speed of the identifying vector field is decreased by a 
change of variable: let r be the function defined by: 



[0,1] — > [0,2] 
* s i ^ -2(Vi-l). (2 ' 3) 

is then given by s = and r = 2 and Eq 0> is given by s = 1 and r = 0. The new identifying 
vector field l/ a is then chosen such that: 

dr(V a ) = 1 so that V a = (r*-R)'(l - 2mR)J^-d R . (2.4) 

y M 

The foliation JC S , when parametrized by r, is denoted by E T . 
Finally, we can prove the following estimates: 

Proposition 2.7. The following equivalence holds: 

E(S+)+E(S U0 )*E(X u °<) 

where 

E(S+)= [ i* f+ {*f a T ab ) and E(S Uo ) = f i* s (*f a T ab ). 

Proof. The proof of these estimates relies on Stokes theorem applied between the hypersufaces 
S^ Q = {(u,R,wgt)\u — uq, \u\ > sr*}, t K s ^ T ) = S T and Eq 0> . This theorem can be used here 
since the data are compactly supported in So and, as a consequence, the future of the support 
of the initial data does not contain the singularity i°. Let be M* o be the subset of fi+ whose 
boundary consists in these hypersufaces. We have: 

/ i* Su (*T a T ab ) + f *l T (*T a T ab ) - / i£ o (*T«T ab )- = [ V" (f b T ab ) M [.g] 

and, using the notations in the proposition, the foliation given by r defined by equation (|2.3p 
and lemma |2~51 this becomes: 

E(S^)+ [ il T {*f a T ab )~ [ z| (*f a T ah ) 
4mi? 2 (3 + uR) {d R (f)) 2 + (1 - 12mR) 4> (u 2 d u (j) - 2(1 + uR)d R <f>) 

+f a V a b—\(r*R)^(l - 2mR)J -^-duAdugi ) dr 
4 J \j \ u \ J 

The error term is bounded above in absolute value; each term is evaluated separetly: 
\{r*R)Hl - 2mR)J^-AmR 2 {2, + uR) (d R (/)) 2 | 



■u 



*R)i(l - 2mR)(3 + \u\R)(R\u\)^ 2 R^- (d R 



u 



<(l + e)l.l.(4 + e).(l + e) ^(^) 2 
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and 



| (1 - VlmR) (r*i?)5(l - 2mR) 

1 



(1 - \2mR) - 2mR)^/R\u\4>{ud u 4>\) 



< -(1 + 6e)(l + e)» • 1 • (1 + e) (0 2 + (ua u </)) 2 ) 

< 2 + ( U a u ^) 2 . 



The remaining term is controlled by: 

/ /? 

(1 - 12mR){r*R)%{l - 2mR)J — 



Remark 2.8. This term is the main obstacle in the use of the parameter s: if the foliation was 
parametrized by s, this term would be replaced by: 



< (l + 6e)(l + e)V| d^d R 4 

u\ 

< (l + 6e)(l + e)* I <p z + ^-Ad_ 



. (r*i?) 2 (l - 2mR) ,„ , 
(1 - 12mRy L <j>d R <j> 



which cannot be compared to <f) 2 + T^r(^fl<^) 2 
Gathering these inequalities, it remains: 



< (l + 6e)(l + e) 2 



4mi? 2 (3 + uR) (d R (f>) 2 + (1 - !2mR) (/> (u 2 d u (j) - 2(1 + uR)d t 



R 



T a V a b— Ur*R)i(l-2mR)J—duAduj S 2 dr 



< 



E(E r )dr. 



Finally, the following inequalities hold: 



£(E r ) + E(S& t) ) < E(S r )dr + E(I%») 
E(Zo 0> ) < Jo E(Zr)dr + E(E T ) + E(S& r) ) 



(2.5) 



Since E(St^) 1S positive, the integral inequality holds: 



E(Z T )< f E(j: r )dr + E(^ 0> ). 
Jo 



Using Gronwall's lemma, we obtain: 

B(S T ) < E(^ 0> ). (2.6) 
Putting this inequality back into (|2.5p . we obtain the first part of the inequality, for r = 2: 

E(S+) + E(S Uo )<E(Z u °>). 

The other inequality is obtained by doing the same calculation between the hypersurfaces 
S Uo , s = {(u,R,lo S 2)\u = u , \u\ < sr*}, J( s and ^ 0> . Let be M s u ° be the subset of fi+ whose 
boundary consists of these hypersurfaces. Applying Stokes theorem, using the parametrization 
by t and the previous estimates of the error term, we obtain: 



JO 



E(Z r )dr + E(J?+)+E(S s V) 
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As a consequence, since the integrand in E(S U °' S ^) is positive, the following integral inequality 
holds: 




□ 



2.2 Energy estimates far from the spacelike infinity i° 

The estimates which are obtained in this section are widely inspired by the work of Hormander 
[3"2] and generalized in [33] to establish the existence of solutions for the characteristic Cauchy 
problem for the wave equation on a curved background. The main tool consists in writing the 
characteristic, or weakly characteristic (that is to say is locally either spacelike or degenerate; 
this is also referred to as achronal) surface as the graph of a function and expressing all the 
relevant quantities in term of this graph. This method was also used by Mason-Nicolas in |34) 
to control how spacelike surfaces converge to null infinity. 
M\fi+ n J + (Eo) is divided in two parts as follows: 

• let E be a spacelike hypersurface in M for the metric g such that E n J? + = S Uo n J? + 
and S is orthogonal to T a . 

• the part of M\f2+ contained between So and E, denoted by V; 

• and finally the future of E, containing i + , U. The subset of its boundary in + is denoted 



This decomposition of the future of So is represented in figure 12.21 

Remark 2.9. The fact that the timelike vector field T a is orthogonal to the spacelike hyper- 
surface is an assumption which is made about this timelike vector field. It is built outside the 
Schwarzschildean part of the manifold as follows: 

• outside fi+ , T a is smoothly extended such that it is orthogonal to a uniformly spacelike 
hypersurface E whose boundary is ,y + D S Uo . 

• The intersection of S UQ and J? + is the Schwarzschildean part and is given {u = u , R = 0} . 
The vector field T° is then orthogonal to this 2-dimensional surface. 

This subsection deals with estimates on U. As previously said, we use here Hormander's 
technique which consists in writing J? as the graph of a function. 

We consider on U the flow associated with eS = . T ~ ~, . Let t be the time function induced 

u gadT<--T d 

by this flow. The vector field eg is completed in an orthonormal basis of TM by choosing an 
orthonormal basis {ef;i = 1,2,3} of the spacelike foliation {E t } induced by t. For the sake of 
clarity, S is denoted St as corresponding to the slice {t = T} (T is chosen to be non zero, in 
order not to introduce confusion with So). 

2.2.1 Geometric description of 

Using the flow ^> t associated with eg, can be identified with S^: 



by A 




(2.7) 
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Figure 2: Future of Sq 



where ip(x) is the time at which the curve t H> ^(x) hits J?^ ■ <?t can then be considered as 
defined by the graph of the function (p : x G £<r 1 — I l p( x )- 

We denote by Viip the derivatives of ip with respect to the vector tangent to £r e? at time 

T. 

Remark 2.10. a. As noticed in the introduction, the spacetimes constructed by Chrusciel- 
Delay and Corvino-Schoen have the specificity that the regularity at J^ + \{i°,i + } can be 
specified arbitrarily. In order to insure that some geometric quantities are defined, we 
assumed that the manifold is C 2 differentiate at i + }. The use of implicit function 

theorem then insures that the function tp has the same regularity. 

b. The function (p is defined on a compact set and as such admits a maximum. This maximum 
is denoted by T max . 

The lapse function associated with this choice of time t is denoted by N. The metric can be 
decomposed as: 

g = N 2 (dt) 2 -/i St 

where h is a Riemannian metric on E t depending on the spacelike leaves of the foliation induced 
by the time function, and N is the (positive) lapse function. Since the time function is built 
from the vector eg, the vector field dt satisfies: 

d t = Nel 

The following lemma describes the geometry of in term of the parametrization: 
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Lemma 2.11. The vector N a defined by 

N a = e a - 

je{i,2,3} 

is normal and tangent to the hypersurface J'^ . 
The set of vectors defined by, for i G {1, 2, 3}, 

tl = NVweZ - el 

are normal to N a and, as such, forms a basis of T.y^ . 

Proof. The fact that 7V a is null directly comes from the fact N a is normal to , which is a 
null surface. The derivatives of cp then satisfy: 

iV 2 £ (V^) 2 = 1 

i=l,2,3 

The tangent plane to J?^ is given by the kernel of the differential of the application 

x i — > ((f(x),x) 

which is given by 

h a e St i— > h a + gij{ip{x), x)V l V h j Ne a . 

St 

It is then clear that the set of vectors 

*? = TVVVeo - e? 
forms a basis of T,y^ . N a is then a linear combination of them: 

N a = Y NV'iptf. 

i=l,2,3 

□ 

As a direct consequence of this lemma, the vector defined by 

r a = e a + Y NV ^ e i 
j'e{i,2,3} 

is null and transverse to Jfjf . 

In order to complete the geometric description of J 2 ^" and facilitate the calculation afterwards, 
we introduce the following objects: 

• using the Geroch- Held- Penrose formalism, the set of two null vectors (r a , N a ) is completed 
by two normalized spacelike vectors (v^wf) tangent to J"^ to form a basis of TM; 

• .y^ is endowed with the 3-form: 

\tj = t\ A t 2 a A tl 

which satisfies: 

r a AiJ M 2 a Ml = {l + N 2 Y, (VV) 2 M.9] 

i=l,2,3 

= M% 

this 3-form will be used as the form of reference to calculate the energy on J^jf. 
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Remark 2.12. The tangent vector n a used in definition \ 1 . 1 2\ for the H l -norm on Jfi + is colinear 
to the vector N a : 

n a 

N a 



9cdT c T d 

As a consequence, the norms associated with these vector fields are equivalent. 

Finally, in order to prepare the estimates, the expression of the 3-form -kT a T ab on the surfaces 
and is given: 

Lemma 2.13. The restrictions of the energy 3-form *T a T ab to T, t , for t given in [T, T max ], and 
ar e given by, respectively: 



i% (*T a T ab )=\\T a \\ (±5>° V »< 



2 4 



and 

(*f a T ab ) = D (V a V Q 0) 2 + «V a </>) 2 + (v%V a q>f) + ^(>^A t* 

Remark 2.14. a. TTie expression which is given for the energy form on J? is consistent 
with the one obtain by Hormander, since it only depends on tangential derivatives to null 
infinity. It is nonetheless not identical: the result of Hormander is similar to a calculation 
made with respect to the Riemannian metric obtained from g and the timelike vector field 

b. The part of the energy form on J? given by (wfV a ij4) 2 + (uf V a<f>) 2 is usually interpreted as 
the norm of the gradient on a 2 sphere, even though the distribution of 2-planes Span(v\, V2) 
is not integrable. 

Proof. Using the basis (e?)j-o,...,3 which is adapted to the foliation, the energy 3-form over E t 
can easily be calculated: 

i* t (*f a T ab ) = (f a V a 4>) *: St (V b 0) + (- l -g cd V c <pvU+ y + i^(*T b ) 

Since the vectors {ef}j = i.2,3 are tangent to the hypersufaces E t , we obtain: 
. i^ t (*T a ) = \\f a \\e^\g] = ||T Q || e i A e 2 A 

• i%(V b <f>) = e b V b cf>(e$Mg}) = (egVb^eJ A e 2 a A e 3 a 
and 



z* St (*TT ab )=\\t*\\ (I^( e ?V^ 

V i=0 



2 4 



To calculate the restriction of the energy form to J^, the vector fields T a and V</> are split 
over the basis ( N a , r a , , v 1 ^ ) : 

llT b ll 

T b = JJ — _ — LL ( N b + T b^ (2.8) 
V c 0V c = iV a V a 0r Q V a (/)-KV a 0) 2 -( W2 a V a 0) 2 . (2.9) 



21 



The only relevant terms in the expressions of i* Jf+ (*Vb</>) and i*^ + (Tf,) are those which are 
transverse to J^, so that: 

and, finally, using equations (|2.8p and (J2T9J) , since the function b vanishes at J*: 

*f a T ab j = ^V a 07V Q V a - l -V c ctN c 4> + A t\ A t\ 

\\T a \ 



((7V a V a 0) 2 + «V Q 0) 2 + KV a 0) 2 ) + 2 ) ti A tl A £. 



□ 



2.2.2 Energy estimates on U 



The techniques used in this section are exactly the same as in the other section: they rely on 
Gronwall lemma and Stokes theorem carefully applied to the 3- form -kT a T a i,. 

The first step of the proof consists in establishing a decay result of the energy on slices 
{t = constant}. 

Proposition 2.15. Let E(T, t ) be the energy on the slice S t : 

E (^) = I ^E«v Q 0) 2 + ^+4^ St 

where /i Ef = el A e 2 A e\. 
Then this energy satisfies: 

and for s and t two real numbers in [T, T max ], such that t > s: 

Remark 2.16. In this section, the calculations are made with respect to E(Y,t) rather than 
fvjh (*T a T ab )- 

Proof. Since T a is a non- vanishing timelike vector field on the compact M, there exists a positive 
constant C such that: 

^<\\f a \\<c 

and, as a consequence of lemma 12. 131 the energy E(Y> t ) is equivalent to J Ej i^ t (*T a T a bj since: 
Remark 2.17. The same result holds for the energy on /y, that it to say that: 



i 



i*^ (*T a T Q6 ) « ^ + (V a V Q 0) 2 + «V a </>) 2 + KV a </») 2 ) + |^ tl A £ A i* 
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We denote by E(J^^) the right-hand side of this equation. The expression of this energy is not 
intrisic, since it depends, on the one hand, on the parametrization of by the function tp 
and, on the other hand, on the choice of a basis. The energy used by Hormander has the same 
property that it depends on the graph and on the chosen coordinates. 

We assume here that t > s and apply Stokes theorem between the surfaces T, t and S s . The 
part of ^ + between the time t and s is denoted J*l and the part of U between E t and S s , C7*. 
The following equality holds: 



& ( *T a T ab 



i% ( *T a T ab 



i|. [*T a T ab 



J (i/(-f^T ab + (1 - iscal § )#f a V a + f a V a ^j M 



As it was noticed in lemma 12.131 the integral over ,f + of the energy 3-form restricted to ,f is 
positive. So the following inequality holds: 

£ i% (*t a T afc ) + i% (*f a T Qb ) - £ (*f a T Q 

< ^ (V a T fc >T a6 + (1 - iscal § )^f a V a + f a V a 6^ 
and, as a consequence, 





V {a f b) T ab 


+ 









(1 - -Scal § )0T a V a ^ 



Since i|; is non- negative (see remark [2. 17[) and 



T a V a b 



it remains: 



/ 


y{afb) Tab 


+ 









i% (*T a T ab ) «£(£ r ), 



(1 - -Scalg)</»T a V a( ; 







+ 


f a V a b 



-li,[g\+E{Y, s ), 



Since U is compact, there exists a contant c depending on V Q T b , Scalg and T a \7 a b which controls 
each term in the error 

J (V a T fc >T ab + (1 - iscal § )0f a V a + f a V a b^j fx[g] 

in function of the energy on a slice at time r. To perform such an estimate, the 2-form \/\ a T b > 

is split over the orthonormal basis (e?)i = o 4; each of the components is bounded by c. The 

remaining terms, when contracting with T ab , are either products of derivatives or functions which 
can be estimated by <fi 2 or b<p 4 and their derivatives. 

Finally, the volume form fi{g] is decomposed over the basis (e")j=o,...,4 to obtain: 



£(£*)< J E(E r )dr + £(£ s ) 
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where the form dr is e a . So, applying Gronwall lemma, we obtain, since we are working in finite 
time: 

E(E t ) < £(£,). 

□ 

A straightforward consequence of this proposition is that all the energies on slices are con- 
trolled by the energy on E-r- This is a necessary step to establish the following proposition: 

Proposition 2.18. The energies on J^f and on Ey o,re equivalent: 

E(y+) a £(£ T ). 

Proof. The proof is based on the use of Stokes theorem. We denote by, for t between T and 
T 

± max ■ 

• Ut the part of U for time greater than t; 

• J^ t + the part of J :+ for time greater than t. 

Stokes theorem is used between the hypersurfaces ^ + and E t : 

[ Ut a T ab ) - / i% Ut a T ab ) 

= (v<- a f b >T ab + (1 - iscalg)0f a V a + f a V a ^pJ fx[g] 
Using exactly the same estimate as in proposition 12.151 of the error term, we obtain: 



< 



E(Z r )dr. 



*T a T ab ) - i% (*T a T ab 
As a consequence, the two following inequalities hold: 

(*T °T ai) ) - £ il t (*T a T ab ) < J £(E r )dr. (2.10) 



and 



since 



£ (*T a T ab ) - (*T Q T ab ) < jf £(E r )dr, (2.11) 



+ (*T°T oi ,J > / (*T a T al 



We first deal with inequality (|2.10p . Since, according to proposition ^. 151 all the energies on 
a slice {t = constant} are controlled by E(St) for t > T, the integral J max E?(E s )dr satisfies: 



£(E r )dr < £(E { ). 
Using inequality (|2.10p . a straightforward consequence is: 

^ + (*f a T ab ) < £(E T ) 

and with remark [2. 171 

£(^ T + )<£(E T ). 
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On the other hand, to obtain the second inequality, we use Gronwall lemma in inequality 
(|2~TTj) ; this gives: 



and, consequently, for t = T: 



2.3 Estimates on V 



E(Z t ) < f Uf a T t 
£(E r ) < E{S+). 



□ 



The geometric situation in this section is almost the same as in the previous one since the 
hypersurface S Uo is known to be null, the only difference being that an additional term comes 
from the boundary of the future of Eq 0< . The energy estimates will then be obtained in exactly 
the same way. Nonetheless, we wish to keep the term with an energy on the hypersurface S Uo , 
in order to compare these terms with the inequalities obtained in section f2. II 

It is clear that the time function which was defined in the previous section cannot be used here 
since T a is not necessarily orthogonal to Eo. We now consider another time function r, defined 
in V (or in a neighborhood of V, as done in remark 12.191 below) such that the hypersurface 
{r = 0} corresponds to Eo and the hypersurface {r = 1} to E = St- We consider on V the 
orthonormal basis (eg, ef)i— 1^2,3 such that: 



ff c dV c TVV 



By construction, the vector fields (e")i=i 2,3 are tangent to the time slices. 
We introduce the following function: 

a = 1 + (9cdf c ef) 2 > 1 

i=l,2,3 

where f c is the normalization with respect to the metric g of the vector field T a . 

Remark 2.19. • Such a time function r can be constructed as follows: since M is globally 
hyperbolic, there exists a time function on M . Let f be a time function. Let tyf be the 
flow associated with r and let Vq be the preimage of E on Eo by the flow. We then obtain 
a diffeomorphism defined by: 

V — > E 
x 1 — > 

where <f>(x) is the time at which the curve f 1— > (x) hits E. The new time function r is 
then defined as: let p be a point lying between Vq and E ; p is written tyf(x); then 

i \ f (p) 

4>{x) 

satisfies the required assumption. 
• T a and eg are both uniformly timelike. Therefore, the scalar product 

(3 - g c J c e a 

defines a positive function over V since T a and eg are both future directed and timelike 
over a compact set. 

The following notations will be used, in order to be coherent with section 12.11 
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• the section of the initial data surface below uq is denoted Eg 0< ; 

• as previously introduced, the slices {f = contant} in V are denoted E t ; 

• the part V between time t and s (with t < s) is denoted by V t s ; 

• the part of S uo between time t and s (with t < s) is denoted by 5|; 

The expression of the energy on S Uo are the same as the one defined in 12. II fsee proposition 12. 31) . 

Following the method already used in this paper, a geometric description of the energy 3- 
form is given and an equivalence result of the integral of the 3-form with an well chosen energy 
is established: 

Proposition 2.20. The restriction of the energy 3-form to E t is given by: 

f (/ C V C 0) 2 



and, as a consequence, the following equivalence holds, for all t in [0, 1] : 

f ll t (*f a T ab )^ [ l(/ C V c 0) 2 + ^ ( e i^a4>) 2 + Y + ) e l A e 2 A e 3- 

•^ St ^ St V i=l,2,3 / 

FFe denote by E(Y, t ) this energy: 

mt)=f ((/ c v c 0) 2 + ( e W) 2 + y + ^I e ? Ae 2 A 4 

•^ E * V i=l,2,3 / 



Remark 2.21. T7ws proposition together with vrovosition \2. 15\ states that the energy on a space- 
like slice for two uniformly timelike (for the metric g ) vector fields are equivalent. This justifies 
that we write in the same way the energy in proposition \2.15\ and in this proposition. 

Proof. The proof is essentially the same as the case where the vector field T a is orthogonal to 
the foliation. It is given in the appendix (see proposition IA.6jl □ 

From now on, the strategy is exactly the same as in subsection 12.2.21 The fact that the 
energy over E dominates the energies on the slices E 4 is established: 

Proposition 2.22. The following estimate holds, for all t in [0, 1]; 

£(E t ) < £(E) 

Proof. Let t be in [0, 1]. The energy 3-form is integrated over the surfaces E, E t and Si; Stokes 
theorem then gives: 

f Ut a T ab ) - f i* Suo Uf a T ab ) - f i% Uf a T ab 
= (v {a f^T ab + (1 - iscal § )#f a V Q </> + f a V a ^pj fi[g 
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As in the proof of proposition 12.181 the error term is estimated considering that V^T 6 ) has 
bounded coefficients in the given basis, that Scalg is bounded and using the behavior of b. The 
volume form is decomposed on the basis (dt, ejji,2,3 as: 



Ms] 



dt A el A e 2 a A e 3 a 



egV a i 



where eg V a i is a positive function and, as a consequence, is bounded above and below by two 
positive constants. 

We obtain, when contracting the volume form with eg in order to integrate in time: 
J vi (y {a f b ^T ab + (1 - iscal § )0T a V a + T a V a ^ fx[g] < E(Z t )dt. 
Using the energy equivalence proved in proposition 12.201 the following inequality then holds: 

E(Sl)+E(Z t ) < J £(E t )di + £(E). (2.12) 

Since E(S}) is non negative (see lemma [L5l and proposition ^. 3|) . f|2 . 1 2[) turns into the integral 
inequality: 



E{^ t )<J^ £(£ 4 )dt + £(£). 



Using Gronwall's lemma, we get: 



£(£ t ) <£(£). 



□ 



Finally, the following proposition holds: 

Proposition 2.23. The following estimates are satisfied on V : 

E(£) * E(S U0 ) + E(Y,%°<) 

Proof. The energy 3-form is integrated over the surfaces Eo, E< and Sq, the application of Stokes 
theorem gives: 



i% (*f a T ab ) - J i* SuQ (*f a T ab ) - [ i* So (*f a T a 



= (v {a f b ^T ab + (l - iscalg)0T a V a ^ + f a V a ^ fx[g}. 

The error term satisfies the same inequality as in proposition 12.221 

(v (a fVT ab + (1 - iscal § )0f a V a + T a V a ^ n\g] < jf E(E s )ds 

As a consequence, the two following inequalities hold: 

£(£*) + £(E ) < f £(£ t )dt + £(£ t ) 
Jo 

E(E t ) < [ E^dt + EiS^+E^o). 
Jo 



(2.13) 
(2.14) 
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The right-hand side of inequality (|2.13[) is estimated via proposition 12.221 as: 

/ £(E 4 )di < £(E) 
Jo 

and, as a consequence, for t = 1, the first part of the equivalence can be stated: 

E(S Uo ) + £(E ) < 
Using the positivity of E(Sq), inequality (|2.14p becomes: 



£(E t )< f E(E t )dt + E(S Uo )+E(Z ). 
Jq 



/o 

Using Gronwall's lemma and setting t = 1, the second part of the equivalence is obtained: 

E(T.) < E(S Uo ) + E(E ). 

□ 

2.4 Final estimates 

Finally, using the three propositions 12.71 12.181 and 12.231 the following a priori global estimates 
hold: 

Theorem 2.24. Let u be a smooth solution with compactly supported data of the nonlinear wave 
equation: 

1 

□it H — ScaLu + bu = 0. 
6 

Then, the a priori estimate holds: 

E(^)^E(J+). 

Remark 2.25. a. The solution of the wave equation is assumed to be smooth in order to 
avoid the problem of defining trace operators for weak solutions of the equation. Nonethe- 
less, using a usual trace theorem, as soon as u is in H^{M), its trace over ^ + and Eo is 
well defined. 

b. Furthermore, in the framework of a (characteristic) Cauchy problem, it is know that the 
solution is in TJ^M). Using the same theorem of existence of trace operators, u is only in 
H2(Eq) or Hz(.y + ), which is not sufficient to write such estimates. It will be shown in 
section [4-l\ that these operators are well defined and with values in ii^Eo) or H 1 (J^ + ). 

Proof. Let us consider the hypersurface Eq. This hypersurface is split in Eq 0< and Eq 0> : 

£?(S ) =£(E " 0< )+i;(E " 0> ). 



E(^>)<E(y+)+E(S U0 ), 
E(S U0 )+E(^°<)<E(Z T ), 



Using proposition 12. 71 
proposition 12.231 
and proposition 12.181 

£(E T ) < E(J?+), 
we obtain the first part of the apriori estimate: 

E(E ) < E{S+) + E(S+) = E(S + ). 
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Conversely, let us consider 
Using proposition 12.181 

E(S+)<E(Er), 

proposition 12.231 

E(X T )<E(S U0 )+E(^<), 

and proposition 12.71 

E(S U0 ) + E(J?+) < E(Y,%°>), 
we get the other side of the inequality: 

E{y + ) < E(Z% 0> ) + E(Z% 0< ) = #(£„). 

□ 

3 Goursat problem on J? + 

We show in this section that there exists a unique solution to the Goursat problem on J? + for 
characteristic data in H 1 (,/ + ): 

□</>+±ScaL j </> + &0 3 = ( 

It is known (see [32]) that the linear Goursat problem on a smooth weakly hypersurface admits 
a global solution; nonetheless, due to technical problems coming from the singularity at i° 
(essentially the fact that some Sobolev embeddings are not valid), the existence of a global 
solution must be justified carefully. 

The proof of the existence for problem (|3.1[) is made in three steps: 

a. considering two solutions of the wave equation, estimates on the difference of these solutions 
are established; the main technical problem, which consists in obtaining uniform Sobolev 
estimates of the nonlinearity, is encountered and solved in section 13.11 The proof is based 
on an estimate of the norm of a certain Sobolev embedding and on the decay assumption 
near V s1 of the function b. 

b. Let § be a uniformly spacelike hypersurface for the metric g in the future of So and close 
enough to + . The existence of solutions for problem (|3.ip with characteristic data whose 
compact support contains neither i nor i + is obtained in the future of § for small data in 
section f3. 21 by slowing down the propagation of waves. 

c. Then this solution is extended down to So, by solving a Cauchy problem on § and using 
a density result in section 13.3.11 

d. Finally, using estimates for the propagator between .y + and So obtained in section 13.11 
the result is extended to H 1 (J? + ) for small data in section 13.3.21 

3.1 Continuity result 

This section is devoted to the proof of a continuity result in function of the characteristic data, 
although an existence theorem has not yet been stated. Consequences of these estimates will be 
required to obtain the well-posedness of problem p.l[) . 
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3.1.1 Sobolev embeddings over a foliation 

One of the main problems when dealing with a characteristic Cauchy problem for a nonlinear 
wave equation is the fact that, when using a spacelike foliation of the interior of a light cone, 
with leaves that shrink when approaching the vertex of the cone, the nonlinearity cannot be 
controlled uniformly on these leaves via Soblev embeddings. More precisely, it is indeed known 
(see |31j ) that, when dealing with a 3-dimensional Riemannian manifold with boundary (M,g), 
the following inequality holds: there exist a constant K, depending only on the dimension of the 
manifold, a real number e, and a constant B, such that, for any function in H 1 (M): 

IMIie (M) <(K + e)\\Vu\\ 2 L2{M) + B\\u\\ 2 L2{M) (3.2) 

where the second best constant B satisfies the inequality: 

B > (Vo\(M,g))~i . 

As a consequence, when considering a foliation whose leaves are shrinking to a point, it is not 
possible to have uniform Sobolev embeddings over such a foliation. 

This kind of embedding was thouroughly studied by Miiller zum Hagen et al. (see for instance 
[TU] and [37]) and extended by Dossa (see [10]) with applications to the characteristic Cauchy 
problem for different kinds of equations (system of scalar nonlinear hyperbolic equations in [22] , 
Einstein- Yang-Mills equations in [21], for instance). M. Dossa especially states that, following 
Miiller zum Hagen, the Sobolev constant associated with the embedding of H 1 into L 6 satisfies 
the following result (theorem 2.1.4 in [20], adapted to dimension 3): 

Proposition 3.1. Let us consider in the Minkowski spacetime R 3,1 the cone defined by 

= \J (x 1 ) 2 + (x 2 ) 2 + (a: 3 ) 2 j-. Then, there exists a constant C such that, for any t > and 

for any function u in H 1 (Y^ t ), where H t is the intersection of the interior of the cone and the 
time slice at time t, the following Sobolev embedding is satisfied: 

C.. 

|M|£6(E t ) < ylMlffi(St) 

The proof of this proposition is obtained by doing an homothety on the Sobolev inequality 
over the unit ball. It must be noted that this inequality is generalized to arbitrary metrics. 

We intend here to prove the same type of inequality for a uniformly spacelike foliation 
transverse to the null conformal infinity ^ + . The first step consists in blowing up the singularity 
in i + and endowing the obtained cylinder with a conformal metric which admits a regular and 
non degenerate extension up to the preimage of i + : 

Proposition 3.2. Let n a be a past oriented null generator of J? + and r be an affine parameter 
associated with n a such that r = in i + . 

Then, for any R > 0, there exists a global diffeomorphism $ : (0, R) x § 2 — > n J + ({r = R})- 
Furthermore, the conformal covering of (j^ + , g\jr+) by ((0,i?) x § 2 , ^ 2 - 1 i>*g\ jz-+) can be extended 
smoothly to ([0,i?) x § 2 , ^t^*g\^+) such that the restrictions of the metric —^2^*g\y+ to {r} x 
§ 2 are uniformly equivalent to the standard Euclidean metric on the two dimensional sphere S 2 . 

Proof. Since the metric on ^ + is the restriction of a smooth (C 2 ) metric defined in an extension 
of M, it is possible to use a result of Hafner-Nicolas in [35] to obtain the result. □ 

Let us then consider a foliation of the interior of by uniformly spacelike hypersurfaces 
(S r ) re r 0j m whose intersections with are the level sets of the affine parameter r. For this 
choice of foliation, the following proposition holds: 
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Proposition 3.3. There exists a constant C such that, for each r £ (0, R] and for any function 
u in iJ 1 (E r ), the following inequality holds: 



< c 



\Wu\\ 2 dfi[g\^} + - 



l dfi\g\i 



Proof. For each r > 0, consider a diffeomorphism <I> r of the three dimensional ball B(Q, 1) into S r 
whose restriction on the boundary is the restriction of the application ^ to the two dimensional 
sphere {r} x § 2 . Such a map exists since the manifold is globally hyperbolic. Using proposition 
13.21 the family of Riemannian metrics — T$*g\-£ r for r £ (0, R] on B(0, 1) is uniformly equivalent 
to the Euclidean metric on the unit ball B(0, 1). 

Let r be in (0, R] and u in H 1 (£ r , -g\s r ). u o $ r is then in H 1 (B(0, 1), -$*g| Sr ) and, a 
fortiori, for a given r > 0, in H 1 (B(0,l),—^Q*g\z r ). Let us consider the Sobolev inequality 
in B(0, 1) for the Euclidean metric: there exists a constant c, depending only on the geometry 
of B(0, 1) such that, for any function / in H 1 (B(0, 1)), the following inequality is satisfied: 



fdx 



< c 



fdx 



>B (0,1) 



'B(O.l) 



l dx 



B(0,1) 



Since the family of Riemannian metrics ~ jz$*g\^ r on B(Q, 1) is uniformly equivalent to the 
Euclidean metric on i?(0, 1), there exists a constant C, which depends only on c and on the 
L°°-bounds of this family of metrics, such that, uniformly in r G (0, R], for any function u in 
H 1 (E r , — <?|i; r ), the following inequality holds: 



(u o $ r ) 6 d/x 



5(0,1) 



< 




(u o $ r ) 2 d/U 



L> r (u o $ r ) || 2 d^ 



where D r is the connection associated with the metric — ^■ < l > *.g|s r . Using the property of the 
conformal transformation and coming back to the manifold S r , it follows: 



u 6 dfi[g\i 



< C 



\Vu\\ 2 dn[g\^ r 



1 



u 2 dn[g\z r ] 



(3.3) 



□ 



Remark 3.4. • It is clear that, using the same arguments as in the previous proof, when 
considering a uniformly spacelike compact foliation whose volume of the leaves has a uni- 
form lower bound, it is possible to obtain in the same way uniform Sobolev estimates. 



It must be noted that formula (|3.3p agrees with equation 
of the Sobolev constants. 



when considering the behavior 



3.1.2 Continuity in term of the initial data for the Cauchy problem 

The purpose of this section is to establish estimates on the difference of two solutions of the 
wave equation. The purpose of these estimates is to obtain continuity in term of initial data, 
characteristic or not. This step is an important one to obtain the continuity of the scattering 
operator. 

Using the method developed below, it is then possible to obtain the apriori estimates of 
section O 
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Let <p an d ip be two solutions of: 

1 , 
Du + -Scahu + bu A = 0. 
6 

We assume that they satisfy one of the problems: 

• an initial value problem on Eo with data in iJg(Eo) x L 2 (Eo) with compact support in Eo; 

• an characteristic initial value problem with data in H 1 (j^ + ) with compact support which 
contains neither i nor i + . 

This ensures that the supports of u and v do not contain the singularity i°. 

Theorem 3.5. Let <p and ip two smooth solutions of the nonlinear problem: 

1 

□u + —ScalgU + bu = 0. 

Then there exist two constants, depending on Scal g , ^j( a T b \ b and the energies of (j) and ip on 
Eo and ^ + such that the following estimates hold: 

\\0-ip\\m { s+) <C(^(E ),^(E ))^_^(E ) 

and 

S _^(E O ) < V'IIhI(jH-), 
where the energy of a function on Eo is chosen to be: 

£ (E O )« f ^ (f a T ab )« f (T a V0) 2 + Y, (e?V Q 0) 2 +0 2 ^ So , 

J Ho JSo i=l,2,3 

where T ao is the energy tensor associated with the linear wave equation and (e^)i = i t 2,3 is an 
orthonormal basis o/TEo- 

Remark 3.6. Because of the a priori estimates, the constants can be chosen indifferently to 
depend on the energy on J? + or Eo. 

Proof. The proof relies on exactly the same strategy as in the first section when establishing the 
a priori estimates. Let 8 be the difference between 4> and ip: 

5 = <f> — f/j. 

5 satisfies the partial differential equation: 

OS + -Scal§<5 + 6(V> 2 + + (f) 2 )5 = 0. 

To establish the inequality, let us consider the energy tensor associated with the linear equation: 

- - f I- - 6 2 

Tab = V a< 5V fc( 5 + g ab --\7 C 5\7 C S 



2 2 



The error term associated with this tensor is: 



V a {f b T ab ) = (V( a f b )T ab ) + 5(f a \7 a 6) - (f a \7 a 5) Qscal^ + b(^ 2 + U + 2 )^ 
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The term A can be estimated by: 

A < I (> + (f "V a <5) 2 ) + \ sup (IScalgl) (<5 2 + (T Q V a <5) 2 ) + 2b ((T a V Q <5)(0 2 + V 2 )*) 

< 2max^sup(|Scal § |),lJ ((f Q V Q (5) 2 + <5 2 + 6 2 <5 2 i/' 4 + & 2 <5 2 4 ) . 

The estimates will then be obtained in exactly the same way as in section 1, provided that 
we are able to use the Sobolev embeddings on the spacelike slices. The main problem then arises 
when working in the Schwarzschildean section because of the choice of the foliation 5f s which 
contains the singularity i°. 

On U, the analogue of the equation (|2.11| is: 



< max 



{l,sup(|Scal § |),sup(||V( a f f ')||)} (£ ma * (EsW+J^ (b 2 5 2 ^ + b 2 6 2 ^ 4 )^dt\ , 



where: 



Espt) = I \ ( E( e i V ^) 2 + y ) ^ = \ (iMlW*) + \\T a V a S\\l^ t 

jTl t \ l=0 J 



(3.4) 



The first problem arises when dealing with the terms: 

S t 



The calculations are made for the first term only. Using Holder inequality, we have: 

fc^W < ( sup b) ( f 6 6 ^ t ) 3 ( [ VVs t 
Using proposition 13. 31 up to a modification of the affme parameter on , it comes: 



s, 



and, consequently, using assumption [d] on function b, the following uniform Sobolev estimate 
holds: there exists a constant c such that: 



^ ^Vs t <c||<5||^ (St) ||^|| 4 ffl(St) . 
The same inequality holds with <j). The integral inequality then becomes: 

< (ci + iHi 4 ffl(s) + imii hs) ) ^ Tmax E 5 (x t )dtj . 

This gives, using Gronwall's lemma: 

E 4 (E t ) < CX P ((d + 11011^^) + HVII^CS)) Pmai - *)) Mm^+y 
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The other estimate is obtained when noticing that the inequality also holds: 

< (Ci + Mfmw + M\H*m) {^[^ Es&t)toj • 

Using proposition 12. 18| the 7J 1 -norm of (j> and ip on £ is controlled by the if 1 -norm of tp and 4> 
on and, as a consequence, by the i? 1 -norm of tp and <j> on ^ + . We then finally obtain on 
U the following inequality: there exist two increasing functions cjj and Cjj such that: 

B*(S) < C [/ ((||0||l 1( ^ +) + ||Vllffi(^ + ))ll^||Hi(^ + ) 

< C £7 ((||^||l rl( ^ +) + ||V||^ 1( ^ +) )^(S) (3.5) 

On V, the principle is exactly the same: the only modification comes from the fact that a 
new boundary term arises, corresponding to the boundary of the Schwarzschildean section. We 
work with the same geometric configuration. The equivalent of equation (|2.12p is then: 

Es(Sl) + £?a(Et) - E S (E) <C 2 I E s {^ t )dt + f b 2 5 2 ip A + b 2 S 2 ^dt 
where E$ has the same expression as in equation (|3.4|) . We finally obtain the energy equivalence: 
£,(£) < c v (f\\<l>\\U ( s +) + {E 5 (Z%°<) + E 5 (S U0 )) 

E 6 (T%><) + E S (S U0 ) < C v (m\mv+) + \\1>\\Hi(s+))E e p) (3.6) 

Remark 3.7. In the subset V of M, the energy on a slice is controlled by the upper slice, which 
is denoted by £ as said in proposition \2.22l As this energy is controlled by proposition \2.18\ 
by the H 1 -norm on ,J?^ and, as a consequence, on , this explains why the energy on 
appears in the inequality. 

Finally, on f2+ o , it is not possible to use the same method as above to control uniformly the 
Sobolev constant. The strategy consists in adopting the same foliation by the hypersurfaces 
JC S . The energy on this foliation is weighted Sobolev norm with a precise decay. The Sobolev 
embeddings must then be adaptated to that decay. The identifying vector field is used to write 
the integral (see formulae (|2.3p and (|2.4p V The error term can then be expressed as: 



r(<0 

o \J3t r 



4mi? 2 (3 + uR) {d R <t>) 2 + (1 - 12mR) <p {u 2 d u <t> - 2(1 + uR)d R q 



-{T a V a 5)5{(f 2 + <pip + if?) }(r*R)t(l-2mR) x l-^ 1 duAdu) § 2 jdr 



In this subset of M, the error is bounded above by, using Holder inequality: 
< [ {S) W) + y=|Z (J 6°du A da*) " (J + du A d. s2 ) ' | dr. 
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< I (d u <t>+ Vi L d R 4>) +|V S 2 0| 2 + 2 dwdw S 2 



< / 2(d u 0) 2 +2(r*R) 2 (l-2mR) 2 ( (d R <f> 2 )\ + V S 20| 2 + tfdudu^ 

J]u , + oo\xS 2 \\ u \/ 



The Sobolev embedding from H 1 into L 6 must then be realized uniformly in r with respect to the 
volume form dwAwg 2 , which is the volume form associated with the cylinder ]uq, +oo[x§ 2 and the 
metric (du) 2 + dw| 2 - Since the Sobolev embedding from H 1 Quo, +oo[xS 2 ) into L 6 (]uo, +oo[x§ 2 ) 
is valid in this geometry, we obtain, in the coordinates (u,cjs?) the following Sobolev inequality: 

f 

3f s J]uq, +oo[xS 2 

< K I (d u (4> )) 2 + |V S 20| 2 + ^dudtJ S 2 

J]«0,+oo[xS 2 X s = {u=-sr*\ 

']tio,+°o[x§ 2 

']«(), + Oo[x§ 

< / 2^(d u 0) 2 + (l + e) 2 -^-^(d R 2 ) + \V S 2q>\ 2 + cj } 2 dudu Js2 . 

■/] M0 ,+oo[xs 2 «o 2m|u | \u\ 

Then there exists a constant K, depending on uq and e such that, uniformly in s: 

Using the fact that (see equation (^.Gp ): 

the following integral inequality holds: 

\E S {M T{S) + E S (S U0 ) - E 5 (X% a> )\ 

<(C + K{E^>) 2 + £ v ,(£^>) 2 )) ( T{S) E 5 (K T )dT 



t(s) 



<(C + K(E^Z ) 2 + ^(S ) 2 )) / E s (K T )dT 

Jo 

and, using the a priori estimates given by theorem 12.241 

\E S (H T(S) ) + E 5 (S U0 ) - E S (T$>>)\ < (c + K (m\m ( ^) + U\\m ( ^))) j^" E 6 (Xr)dT. 
Finally, there exist two increasing functions Cq+ and C fi + such that: 

E S {J+) + E 5 (S UQ ) < c u + o (iHll rl(J r +) + ||VI|^(> + ))^(^ 0> ) 

E S (K 0> ) < C n + o (U\\% H y +) + \\n%^ +) ){MK)+ E ^u )). (3.7) 

Eventually, combining inequalities (|3.5p , (|3.6p and (|3.7p as in section [5] for the proof of 
theorem 12.241 we get the existence of two increasing functions c and C such that: 

£ 5 (£o) < c((||^||l 1( ., +) + |Hliri(^ + ))||5||iri(> + ) 

< G((ll0ll^(^ + ) + llV'll^ ( ^ + ))^(So). 

Because of the a priori estimates given by theorem 12.241 the iJ 1 -norm of <j> and ip on J> + can 
be replaced by the energy of (f> and ip on Sq. □ 
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This result is equivalent to the result obtained by Hormander at the beginning of his paper. 

Finally, as already mentioned, a by-product of this result is the continuity result for the 
Cauchy problem for the nonlinear wave equation on a uniformly spacelike hypersurface S, trans- 
verse to J !+ . The problems are exactly the same: obtaining uniform Sobolev estimates near i° 
and in the equivalent of the region V. The techniques to solve this problem are then exactly 
the same. We are working with functions <j) and ip which satisfy the same assumptions as for 
theorem 13.51 

Proposition 3.8. Let § be a uniformly spacelike hypersurface transverse to . Let <f> and ip 
be two smooth solutions of the nonlinear problem: 

1 , 
□u H — ScaLu + bu = 0. 
6 

Then there exists two constants, depending on Scalg, V a T b , b and the energy of <f> and ip on Eo 
and S such that the following estimates hold: 

S^_^,(§) < C(E < p(T, ), E^(T, ))E^-^(T, ) 

and 

Ej-^CEo) < C{\\(j}\\ H i ( y +) , ||^|| H i ( ^+ ) )£ , 0_ v ,(§), 
where the energy of a function u on a uniformly spacelike hypersurface E is chosen to be: 




'0 ■ 



where T a b is the energy tensor associated with the linear wave equation and (e")i = i.2,3 is an 
orthonormal basis o/TS. 

3.2 Solution of the Goursat problem near 

This section presents a proof of the local existence of a solution of the characteristic Cauchy 
problem on J? + , which consists in slowing down the propagation speed of waves, as done in 
Hormander 's work in |32| . This gives a weak convergence result. The weak solution is constructed 
down to a certain uniformly spacelike hypersurface, constructed as follows. 

Let us consider a smooth function 9 defined on whose compact support contains neither 
i° nor i + . We consider a uniformly spacelike hypersurface § such that: 

• S is in the future of Eo an d does not contain i°; 

• § is transverse to J !+ \ 

• S is in the past of the support of 9. 

A local existence result for the characteristic data 6 down to the hypersurface § can then be 
stated: 

Proposition 3.9. Let us consider the nonlinear characteristic Cauchy problem on : 

f Du + \ScalgU + bu 3 = 

\ u\^ + = e^H\j+). 

where 9 is a smooth function whose compact support does not contain i + or i . Then this problem 
admits a weak solution in the future of § in H 1 (J + (S)). 
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Proof. Starting from the spacelike hypersurface § whose future contains the support of 9, let t 
be a time function on M with t = on §. The metric is split with respect to this time function 
as follows: 

g = N 2 (dt) 2 -h t , 

where N is the lapse function and h t is a Riemannian metric on the spacelike slice t = constant. 
A real parameter A in (|, 1) is introduced to obtain a new family of metrics gx defined by: 

gx = X 2 N 2 (dt) 2 - ht. 

For a given value of A, the hypersurface is now uniformly spacelike and the initial value 
problem on for the wave operator associated with gx, denoted Q\, is then a Cauchy problem 
on Furthermore, since the support of the data 9 does not contain the singularity i + (and 
the singularity i ), the problem can be considered to be set on a smooth spacelike hypersurface. 
The result by Cagnac and Choquet-Bruhat in [B] can then be used. Let us then consider the 
family of solutions (u\) in H 1 (J + (S)) of the Cauchy problem: 

□ a ua + ^Scalg^itA + bu\ = 
u x = 9 on 
d t ux = on J? + 

The first step of the proof consists in proving that the family (ux) is bounded in H 1 (J + (§)) . 
This can be done by noting that the same inequality as in propositions 12.151 and 12.181 holds for 
the metric gx with £<r = S in this case: as a consequence, there exists a constant cx depending 
on the time function t and depending continuously on the scalar curvature of gx such that: 

E Ux (B)<cx[ i^+(*T a T ab ) (3.8) 

where: 

• the integral on the right-hand side only depends on the function 9; 

• the constant cx is bounded above by a certain constant C since it depends continuously 
on A in [^,1]; 

• and the energy E Ux {§) controls the f/^-norm of Ux- 

As a consequence the family (u\) is bounded in H 1 {J + '(§)). 

Then there exists a sequence (A„) converging towards 1 and a function u in H 1 {J + '(§)) such 
that the sequence (u n = ux n ) converges towards u weakly in H 1 (J + (§)). 

Consequently, there exists a subsequence of (u n ), still denoted (u n ), converging towards u 
strongly in L 2 {J + (2>)). This sequence converges a fortiori strongly in L 1 (J + (§)), using the 
Cauchy-Schwarz inequality, since we are working on a compact neighborhood of i + . 

The next step consists in considering the nonlinear term of the equation. We prove that (m^) 
is a Cauchy sequence in L 1 (J + (§)). Noticing that, for any (n,p) in N 2 

Il3 3ll ^ f \ \ Il2 ii ||2 All II 

\\ u n ~ u p\\H(J+(S)) < 2 ^II u ™IIl 4 (J+(S)) + II u pIIl 4 (J+(S)) J \\ u n - u p I Ii 2 ( J+ (S)) 

and, using the Sobolev embedding from H 1 into L 4 (in dimension 4), there exists a constant C 
such that: 

\\ u l - u plUi(./+(S)) < Cswp [\\Un\\ 2 H irj+( S )\) \\u n ~ U p \ |l 2 (,7+(S))- 

As a consequence, (u^) is converging strongly in L 1 (J + (§)) towards u. 
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It finally remains to prove that the function u satisfies the given characteristic Cauchy prob- 
lem. Consider the trace of the functions (u n ) on JP + . Since (u n ) converges weakly in H 1 (J + (§)) 
towards it, using theorem 14.11 the restriction of (u n ) to converges weakly in H?(J !+ ) to- 
wards the restrictions of u to . Using Rellich-Kondrachov theorem, passing to a subsequence 
if necessary, the constant sequence (u n \y+ = 9) converges towards u\y+ strongly in L 2 (,f + ). 
The function u then satisfies the initial condition u = on J r+ . □ 

3.3 Global characteristic Cauchy problem 

A global Cauchy problem is finally derived in two steps: 

a. a preliminary result about the Cauchy problem for a hypersurface in the future of Eo 
whose past contains i ; 

b. the characteristic Cauchy problem is then solved for small data with compact support 
which contains neither i° nor i + and then extended to functions in H 1 (J f+ ). 

3.3.1 Global Cauchy problem for compactly supported data 

Starting from the same data 6 in H 1 (J^ + ) whose support contains neither i + nor i°, the solution 
obtained in theorem 13. II is extended to the future of So by the means of density results and con- 
tinuity of the propagator. The purpose of this section is to show that is it possible, starting from 
the hypersurface § in the future of Eo, to obtain a solution down to Eo despite the singularity 
in i°. 

Proposition 3.10. Let V a be an orthogonal and normalized vector field to the uniformly space- 
like hypersurface §. 
The non-linear problem on S: 

Dv + \ScalgV + bv 3 = 
< t;| =ee J ff 1 (S) 

admits a global unique solution down to Eo in C° (M., Hq (§)) . 

Proof. The method consists in approximating the solution by solutions of the same problem 
with truncated data since the existence result for the Cauchy problem given by theorem 11.171 
cannot be applied here directly because of the singularity in i°. The uniqueness directly conies 
from theorem 13.51 and its corollary. 

Let (xn)neN be a sequence of smooth functions with compact support in the interior of § 
such that: 

\fn e N, supp (xn) C supp (xn+l) and (J supp (x„) = S\<9S. 

Let (un)neN be the sequence defined by: 

Uv n + iScal§^„ + bvl = 
«»| 5 = Xn£eH 1 (S). 

V a V a v n \ s =XnCtL 2 (S). 

Since the data are with compact support in the interior of S, their pasts do not intersect J? + 
and, as a consequence, i . 

Using proposition 13.81 this sequence converges towards a function v in the past of § down to 
Eo for the L 00 !! 1 norm. This function satisfies the initial conditions: 

v\ s = £ on § and U Q V a w| s = C- 
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Furthermore, proposition 13.81 also gives convergence in _fT 1 (J _ (S)) and, as a consequence, us- 
ing Sobolev embeddings, in L e (J~ (§)). v then satisfies the nonlinear wave equation in the 
distribution sense. □ 



Remark 3.11. A direct consequence of this construction is that the trace of the solution of the 
Cauchy problem on So is in Hq(^q). 

3.3.2 Global characteristic Cauchy problem for initial data in H 1 (,/ + ) 

A global solution to the Goursat problem with compact support which contains neither i + nor 
i° is then obtained by gluing solution of the local characteristic Cauchy problem obtained in 
proposition 13.91 and the solution of a well-chosen Cauchy problem on S: 

Proposition 3.12. Let u be a solution to the Goursat problem for data 9 with compact support 
which contains neither i + nor i . 

Then u can be extended from S down to So in C°(R, iJ 1 (§)). 
Proof. Consider the Cauchy problem on S: 

Dv + iScalgU + few 3 = 
u|s = u e H\§). 
V a V a v\ s =V a i7 a ue L 2 (§). 

According to proposition l3.10l this problem admits a global solution v down to So in C°(R, iJ 1 (S)). 
Finally, the function w defined piecewise by: 

w = u on J + (S) and w = v on J + (S ) n J~{S). 

satisfies the Goursat problem: 

( Ou + iScaLjU + bu 3 = 
I <l>\s + =9EH 1 (<S+) 

□ 

Using proposition 13.121 and the continuity result, we can state the theorem of existence of 
the Goursat problem: 

Theorem 3.13. Let us consider the nonlinear characteristic Cauchy problem on J^ + : 

( Du+ ^ScalgU + bu 3 = 

Then, this problem admits a global unique solution down to the future of So in C (M, H (So)). 

Proof. The proof relies on the density of data with compact support in H l (j? + ) which contains 
neither i° nor i + in H 1 (J^ + ) (proposition II .13)1 and proposition ^. 121 □ 

Remark 3.14. As noticed above, the trace of the solution of the Goursat problem on So is in 
#o(£o). 

4 Construction of the scattering operator 

The construction of the scattering operator can now be done by solving Cauchy problems on 
and Sq via the composition of trace operators. 
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4.1 Existence and continuity of trace operators 

The purpose of this section is to define trace operators for the solutions of the wave equation on 
the hypersurfaces Eo and J?^ . A similar construction can of course be realized on the past null 
infinity . 

These trace operators are obtained using the following theorem ([12], p. 287): 

Theorem 4.1. Let M be a smooth compact manifold with piecewise C 1 boundary and consider 
the application T defined by: 



T 



C°(M) — -> C°(dM) 



f — > / 



dAI ' 



Then, for all s > i the operator T extends uniquely to a continuous map from H S (M) into 
H s -?(dM). 

Existence theorems 11.171 and 13.131 give solutions to the initial (characteristic) problem in 
H 1 (J + (So)). As a consequence, their traces on So and J? + are respectively in H^(Y,q) and 
Hz(J? + ). Nonetheless, using the a priori estimates, they are in fact H 1 {^ + ) and H 1 (T,q). 

Remark 4.2. The singularity in i + is not a threat to the existence of a trace since the manifold 
and the metric can be extended with arbitrary regularity in a neighborhood of i + . The prob- 
lem with the singularity i is avoided since the function spaces H 1 (j^ + ) and Hq(Y.q) are the 
completions of smooth functions whose compact support does not contain i° . 



Let us consider the trace operators: 

C °°(S]o) x qf(£ ) — > HHy+) 

where (f> is the unique solution of the problem: 



T + ._ i k^oj y^oj 7 " \^ j ^ ^ 



D(f>+ iScalg(/) + 60 3 = 
<4 b = 9 G C °°(S ) 

f n V Q 0| = 0eC§°(£ o ) 



obtained by theorem 11.171 and 



T o _ / * ~ ► J ff 1 (So)xL 2 (E ) 

I 9 (^l So .(^v a 0)| So ) ^ 

where $ is the set of smooth functions with compact support which contains neither i + nor i° 
and 4> is the unique solution of the problem: 

U(j) + ±Scab> + bcf = 

obtained by theorem 13. 131 

These operators can be extended to ifg(Eo) and H 1 (J? + ): 

Proposition 4.3. The operator T^ can be extended to a locally Lipschitz operator from Hq(T,q)x 
i 2 (E ) to H 1 (J+). 

The operator can be extended to a locally Lipschitz operator from H l {^ + ) to Hq(Tiq) x 
L 2 (£ ). 
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Proof. The proof is done for the operator Tq (it is exactly the same on the other side). 
Using theorem 13. 51 this operator satisfies: 

v(e,e,te) e (Q°(£ )) 4 , 
\\T+(e,8) - T +(£,£)|| ffl( ^ +) < c(R) (\\e - e||| 1(So) + \\e - f||£ a(Eo) ) 3 . 

As a consequence, since the smooth functions with compact support in So are dense in HqCEq), 
it admits a unique locally Lipschitz extension from i?o(£o) x L 2 (Yjq) into H 1 (J? + ). 

The same proof holds for the operator T?. □ 

As already noted at the beginning of this section, a similar construction can be achieved on 
the past null infinity: there exist two Lipschitz trace operators T_ and T _ defined by: 

f H 1 (S )x L 2 (E ) ff 1 ^) (43) 



where </> is the unique solution of the problem: 

D(j)+ iScalg0 + 6(/> 3 = 

^ f a V a 0| So = 0el 2 (E o ) 

obtained by theorem 11.171 and 



, H\y~) i/ 1 (S )xL 2 (So) 



where </> is the unique solution of the problem: 



U(j) + iScal^ + bcj) 3 = 



obtained by theorem 13.131 



4.2 Conformal scattering operator 

Finally, the conformal scattering operator is obtained as the composition of two trace operators. 
Following the idea of of Friedlander in [23] and applied by Mason-Nicolas in [33] for the Dirac 
and wave equations, the conformal scattering operator S is defined by composing the operators 
T° and T+: 

H\y-) (4.5) 

>H\jr-) (4.6) 
Finally, the following existence result for the conformal scattering operator can be stated: 

Theorem 4.4 (Scattering operator). The operator S is an inversible, locally Lipschitz operator 
from if 1 (=# _ ) into in H 1 (J ;+ ). This operator is called conformal scattering operator. 

Proof. The proof is an immediate consequence of proposition 14.31 □ 

Remark 4.5. The conformal scattering operator was introduced to avoid the use of the spectral 
theory which requires the metric to be static. It is nonetheless possible to talk about geometric 
scattering at least in the Schwarzschildean part of the manifold and wonder wether it is possible 
to establish an equivalence in this region. Some answers to this question can be found in \31$ 
(section 4-2) for the Dirac and Maxwell equations. 



S = T+oT»: H 1 (J 



"0 . Tjll 

and its inverse is given by 

S- 1 = T - oTl : H 1 (y + ) 
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Concluding remarks 

There exist several possible extensions to this work: 

• the case where the metric in the neighborhood of i is the Kerr-Newman metric; 

• the nonlinearity could be modified and the equation could, for instance, be quasilinear, or 
satisfy the null condition; 

• following [35] , these results could be extended to peeling results for the same cubic defo- 
cusing wave equation. 

One of the main problems of general relativity is the construction of solutions to Einstein 
equations. One intermediate step could be to establish the same kind of result for the Yang-Mills 
equations. 



A Appendix 

As previouly mentioned, this appendix is divided in three sections. The first one contains the 
proof of lemma 11.141 The second section gives the necessary details leading to the a priori 
estimates. The last one presents an alternate proof, for small data, of the well posedness of the 
Cauchy problem on conformal infinity. 

A.l Proof of lemma 11.141 



Lemma A.l. Let us consider the set of smooth functions defined in B(0, 1) C K. with support 
which does not contain 0. Then this set is dense in H 1 (B(0, 1)). 

Proof. It is sufficient to prove that the constant function 1 can be approximated by a sequence 
of smooth functions whose compact support does not contain 0. 
Let / be the function defined on R + by: 

• / is a smooth function on R + with value in [0, 1]; 

• / = 1 in [|»+oo); 

• / vanishes in [0, |]. 

Let us consider the sequence of smooth spherically symmetric functions defined by: 

Vn e N,Vx e B(0, l)Mx) = f{n\\x\\). 

They satisfy, for all n in N: 

. ^ n = 1 in B(0, 1)\B(0, £); 

• ipn vanishes in B(0, ^); 

• ip n is a smooth function on B(0, 1) with value in [0, 1] since it vanishes in a neighborhood 
of zero. 

Finally, the difference (1 — ip n ) n converges towards in H ^norm: 

|l-Vn|&i = / ((1- f(nr)) 2 + n 2 (f (nr)) 2 ) ^drd^ 

l 




(1- f(nr)) 2 dr + n 2 J (/' \nr)) 2 r 2 dr 



i] 
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The remaining integral converges towards by Lebesgue theorem. As a consequence, the se- 
quence (ip„) n converges towards the constant 1 in H 1 (B(0, 1)). 

Finally, let / be a function H 1 (B(0, 1)) (or in Hq (B(0, 1))). Cauchy-Schwarz inequality 
gives, for all n in N: 



11/(1-^)11^ 



-^«)ili» + ||/iv(i-^)i||^ + ||iv/i(i-^ n )||^ 



< 2 



Nll-Vn" 2 



H 



f/ 1 



(fi/j n )n is then a sequence of functions in H 1 (B(0, 1)) whose support does not contain which 
converges in H 1 towards /. □ 

A. 2 A priori estimates — Some details 

Most of the calculations which are made in this paper was already made in the papers of Mason- 
Nicolas [Ml I3"5] in the linear case. Nonetheless, for the sake of self-consistency, it has been 
decided to give more details about the proofs of these inequalities. Two kinds of details are 
given here: the first one is about the geometry of the Schwarschild spacetime and the second 
one about the estimates themselves. 

Proposition A. 2. The energy 3-form, written in the coordinates (R,u,9,ip), is given by: 



*T a T ab 



\d u (p) 2 + R 2 (l - 2mR) (u 2 d R <pd u <p - (1 + uR){d R <p) 2 ) 



1 



<P 2 A A 



l - ((2 + uRf - 2mi?V) (c^ 2 ) + u 2 Q|V S 



sin(0)efu A d9 A dip 

j2 jA 



2 , * n 4 4 ^ 

- d> 2 <£ 4 \1 
+ Y + b 4~) sin ( ) dRA d6A d ^ 

+ sin(0) [u 2 d u (p - 2(1 + uR)d R <£\ {-d e (pdu A dR A dip + d^cpdu AffiA d6) 
The restriction of the energy 3-form can be written: 
• to "K R : 



i* Ks {*f a T ab ) = (u 2 {d u <P) 2 +R 2 (l- 2mR)u 2 d R( pd u (l> 

x /(2 + uR) 2 mu 2 R 3 , ,\ ,„ 

+R 2 {1 - 2mR) h ^—t (1 + uR)\ {d R <p 2 ) 

u 2 R 2 {l~2mR) + 2(i + \ Q |Vg20| + + b ^\\ s - m{e)du Ad9A #; 



• to S u : 



i* Su (*T a T ab ) = I - ((2+ uR) 2 - 2mi?V) (d R 4> 2 ) 



v> -IVs 



1 2 4 



sin(0)di? A d0 A dip; 



z> + (*T a T afc ) = u 2 ^) 2 + || v s 



+ 4> z + 6y ] sin(6»)dM A d9 A dip. 



4.°, 



Proof. For the calculation to come, let us denote by A, the quantity: 

1 - « <b 2 (b 4 

We calculate the general form of -kT a T ab : 

*f a T ab = *(u 2 dl - 2(1 + uR)d a R )T ab 
= (u 2 d u (f> - 2(1 + uR)d R (j>)-kV b (j) + u 2 Ad c ^n[g] - 2A(1 + uR)d R ^[g\. 

Using the expression of the volume form in the coordinates {u, R,w§p.) (equation (|2.2p ~). we 
obtain: 

d u j^[g] = di?Ad 2 w S 2, dgj/j,[g] = sin(#)di? A du A d9, 
d R Afj,[g] — -duAd 2 u) S 2, d^A^g] = — sin(#)di? A du A d<p. 

V^> is written in the coordinates (R,u,ujg2) as follows: 

V</> = -d R 4>d u - (d u + R 2 (1- 2mR)d R $)d R - V s ^; 

its norm is then: 

V c </)V c = (-d R cb) 2 g(d u ,d u ) + 2d R ^{d u 4> + i? 2 (l - 2mR)d R 4>))g{d u ,d R ) - |V S 20| 2 
= -R 2 (l - 2mR)(d R dp) 2 - 2d R 4>d u <$> ~ |V§20| 2 . 

The Hodge dual of V a <t> is calculated by splitting V a u over (du, di?, dwg2): 

-kVb4> = * (d u <f>du + dR(/)dR + de4>d9 + d^(f)dil>) 
= d u cbdR A da;§2 — d R 4>du A dw§2 

The gradient of each coordinates is calculated: 

V b u = g ab V a u = -d R V b 6> = -sm(0)d e 

V b R = -d u -R 2 (l-2mR)d R VV = ~d^, 

and, as a consequence, 

*du = du A dw§2 -kd8 = — sin Odu A dR A dijj 

*di? = -dR A dw S 2 + i? 2 (l - 2mR)du A dw S 2 *dip = du A dR A d# 

so that *V a </> is: 

*V a </> = 9 M 0du A dw§2 + (-di? A dw§2 + i? 2 (l - 2mi?)du A dw S 2) 
-<9<9^ sin #du A dR A df/> + S^du A dR A d# 
= (d u (p + R 2 (1- 2mR)d R 4>) du A dw§2 - <9 fl <?(>di? A dw§2 
-dgcj> sin #du A di? A d^ + <9,>du A di? A dd 

The energy 3-form is then: 



*T a T ab 



u 2 {d u <f)) 2 +R 2 (1- 2mR) (u 2 d R (/)d u (f) - (1 + uR)(d R c, 
1, 



+2(1 + uR) ( ^|V S 20| 2 + ^- + ^ 



du A dw§ 2 



l - ((2 + ui?) 2 - 2mi?V) (d^ 2 ) + u 2 Q|V S 2 0| 2 + ^ + 6^ 



di? A dw§2 

+ sin(0) [u 2 d u 4> - 2(1 + uR)d R (p] (~dg(j)du A di? A di/' + d^du A di? A d0) 
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Remark A. 3. To calculate the restrictions of-kT a T ab to the hyper surfaces 5f S; S u and ^ , it 
is necessary to give the restrictions of each of the differentials of the coordinates. Nonetheless, 
since dg and are tangent to 5f s , S u and , the only remaining 3-forms to consider when 
restricting to these hyper surf aces are duA du!§2 and dRA du)§2. This means that only du and dR 
should be taken care of. 



Noticing that 



dr 



dR R 2 {l-2mRY 



we get, on JC S , defined in f2+ by u = —sr* 



dR\ 



R 2 (l - 2mR) 



'K, 



Mi 



t R 2 (l-2mR) 



du\ 



The restriction to JC S of the energy 3-form is then: 

i$t,(*f a T ab ) = (\u 2 {d u ^) 2 + R 2 (l - 2mR) (u 2 d R <bd u cj) - (1 + uR)(d f 



-2(1 + uR) 



R 2 {1- 2mR) 



^({2 + uR) 2 -2mR i u 2 ) (d R cf> 2 ) 



5 T 



du A dcj§2 

i^ s (*f a T ab ) = (u 2 {d u $) 2 + R 2 (l - 2mR)u 2 d R ^d u (j ) 
(2 + uR) 2 mu 2 R 3 



+ 



+R 2 (1 - 2mR) 
u 2 R 2 {l - 2mR) 



2s 



- (1 + uR) (d R ^) 



1 



+ 2(1 + uR)) ( -|V S 2</>| + y + ^-j ) )duAdw S 2. 



is defined by R = 0; so the restriction to + of the energy 3-form is: 



i^+{*T a T ab ) = [ u 2 {d u <PY + [ |V S 2 0| +<j) 2 + ) ) duAdw*". 



Finally, for S u which is defined by {u — constant} , we obtain: 



i*sS*T a T ab 



((2 + uR) 2 - 2mR 3 u 2 ) (d R <j) 2 ) + ■ 



T + 6 T 



dR A dws2 . 



□ 



Proposition A. 4. There exists uq, such that the following energy estimates holds on "K s in 

R 



^0 



:k, 



r (^) 2 + |V S 2< 



— + &— 
2 4 



du A du>g2 



Proof. Let us consider the expression of i^- (-kT a T ab ) 



i*x s (*T a T ab 



u 2 (d u tf 



-R 2 u 2 {l - 2mR)d R (j)d u c l 



-R 2 (l - 2mR) 



[2 + uRf 



2 i? 3 



2s 



(i + uR) (d R <t> 2 ) 



(A.l) 
(A.2) 
(A.3) 



u 2 R 2 (l - 2mR) 



2(1 + uR) 



^|V S 2 ( 



y + 6y ) )duAda; S 2. (A.4) 
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Each of these terms is estimated separately using lemma 11.51 and the obtained estimates are 
summed. Let e be a given positive number and let uo be the non positive constant associated 
to e via lemma [T31 e will be chosen during the proof. 
Nothing needs to be done for (| A. 1|) . 

For (|A.4[) , since u is non-positive and s = — -p- = ^- , we have, on one side: 

u 2 R 2 {l-2mR) +2(i _ ^ A = ^^^^^ _ 2mR ) + 2 (1 - \ U \R) 

< (Rr*){R\u\)(l-2mR) + 2 

< (l-e)(l + e) + 2 

and, on the other side: 

'u 2 R 2 (l - 2mR) 



2(1- \u\R)\ = {Rr*)(l-2mR){R\u\) + 2(1 -\u\R) 

> 1 ■ (1 -e)-{R\u\) + 2(1 -\u\R) 

> 2-(l + e)(J2|u|) 

> 2-(l + e)(l + e) 

> l-2e-e 2 . (A.5) 

e is chosen such as (|A.5|) is positive. 

For (|A.3|> . the proof is slightly more complicated. We have, since u is non-positive and 
s = — tt = on one hand: 

9/ ,/(2 + uR) 2 mu 2 R 3 , \ .„ , n9 

R 2 {\ - 2mR) ( V - ' (1 + uR)j {d R 4>) 2 

= R 2 {1- 2mR) ( r *( 2 ~ H^) 2 _ ( m R)(R\u\)(Rr*) - (1 + uR)] (d R cj)) 2 
\ 2 \ u \ / 

^(« 2 ) (1 - 2mR)(Rr*) (^&- - {m R){R\u\) 2 - M + Mfl H ) (A>6) 

<fo(fl^) a Vi + f^+f!l! + 1 . (1+e; 



Starting from equation (|A.6|) . the lower bound for (|A.3[1 is obtained as follows, setting X = \u\R: 
R (d R cj>f) (1 - 2mR)(Rr*) ( fcM^l! - ( m R)(R\u\) 2 - M + M^j 



u| y \ 2 r* r* 

> ( A(^) 2 ") (1 - e ) ■ ( {2 \ X) \ Rr*) - (mR){Rr*)X 2 - X + X 2 



> (^(^0) 2 ) ^y^ (4 - 6X + 3X 2 6(1 + £)X 2 ) 
The polynomial 4 — 6X + 3X 2 reaches its minimum for X = 1 and equals 1 at J = 1 so that: 
r-(cW») 2 J (1 - 2mR)(Rr*) l± LU_ - (™i?)(i?M) 2 - U + U fl | u 



To deal with (IA.2I). we write 



|i?V(l - 2mR)d R <j ) d u 4>\ = (1 - 2mR) \^R 2 \u\^d R <j\ U^ud u A , 



4G 



so that: 



|i?V(l - 2mR)d R <t>d u <t>\ < (1 - 2mR)\ (I (ud u <f>) 2 + J(ij3| u |3)-^(fc 

z \2 ,3 |u| 

< | W f + i(l + e)3^(^ 
Finally, the following equivalence estimates hold: 
c e J (u\d u 4>f + ^j(c^) 2 + |V S2 0| 2 + ^ + 6^ du A dw S 2 < J (*T a T ab ) 

J i*x 3 (*T a T Qb ) < C £ ^ ^ 2 (9 U 0) 2 + ^ (d R cb) 2 + |V S2 0| 2 + ^ + 6^ du A dco S 2 
where 

C e = max (l, 1(1 + e) 3 , (1 + e) (^-^ + (1 + e)) , (1 - e)(l + e) + 2 



and 



and 



4' 6 W ' 



where P(e) is a polynomial in e. e is chosen such that the constant c e is positive. Using lemma 
11.51 there exists uq, negative, \uq\ large enough, such that the estimates of the coordinates hold 
in and, consequently, the equivalence is true on this domain. □ 

Lemma A. 5. The error is given by: 



V a [T b T ab j = 4mi? 2 (3 + uR) (d B ^) + C 1 ~ l2mR) <i> {u 2 d u <j) - 2(1 + uR)d R <p) + T a V a &-^-. 
Proof. The Killing form of the vector T a is calculated via the Lie derivative of the metric: 

V( a T b ) = Lp<7 

= L f (i? 2 (l - 2mi?))(dw) 2 + 2i? 2 (l - 2mR)L f (du)du - 2L f {du)dR - 2L f (dR)du - 2L f (dw S 2) 
= -4R(1 + uR)(l - 3mR)(du) 2 + 2R 2 (1 - 2mR)d (L f (u)) du - 2d (L f (u)) dR - 2d (L f (R)) du - 
= -4R(1 + uR)(l - 3m,R){du) 2 + 2R 2 (l - 2mR)d(u 2 )du - 2d (u 2 ) dR - 2d (-2(1 + uR)) du 
= -4R(1 + uR)(l - 3mR)(du) 2 + 4uR 2 (l - 2mR)(du) 2 - 4udwdi? + 4(Rdu + udR)du 
= (-2R(1 + uR)(l - 3mR) + 4uR 2 (l - 2mR) + 4R)(du) 2 
= {12mR 2 +4muR 3 )(du 2 ), 

or, in the vector form: 

V (a f b *> = 4mR 2 {3 + uR)d R d R . 
A direct consequence of the above formula is that the Killling form is trace-free: 

V a f a = AmR 2 (3 + uR)g(d R ,d R ) 
= 0. 

The scalar curvature of the rescaled metric is given by equation (|1.2p : choosing O = R gives: 



Seals = R 3 V b V b R 



6 

= 2mR. 
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Finally, the error term is given by: 

v (afb) Tab = 4mR 2( 3 + u E)^ R< f ) y + (l-l2mR){f a W a (l))<f> + T a W a b^ 



4mR 2 (3 + uR) {d R <t>) 2 + (1 - 12mR) <t> (u 2 d u <t> - 2(1 + uR)d R c, 



-T a V a b^j. 



Proposition A. 6. The restriction of the energy 3-form to S t is given by: 

(/ C V C 0) 2 



□ 



i|, (*f a T ab ) 



2(1 + E 4 =i, 2 , 3 (^/ C e d c i ) 2 ) 



id,3 V 1 + ^i=i,a,3 K9cdf c e d CiY ) J 2 4 J 

2 consequence, the following equivalence holds, for all t in [0, 1] : 

/ ij t (*f Q T ab ) » f I (/ C V C 0) 2 + ^ (^V a 0) 2 + £ + 6^ ) e? A e§ A e£ 



and, as a consequence, the following equivalence holds, for all t in [0, 1] 

| (/ C V C 0) 2 + (e?V a 0) 2 + f 

We denote by E(T, t ) this energy: 

Ept)= [ ((/ C V C 0) 2 + i^^) 2 + Y + b T K^Ae^. 

•^ E * y 1=1,2,3 J 

Remark A. 7. This proposition together with vrovosition \2. 1 5\ states that the energy on a space- 
like slice for two uniformly timelike (for the metric g ) vector fields are equivalent. This justifies 
that we write in the same way the energy in proposition \2.15\ and in this proposition. 

Proof. The proof is essentially the same as the case where the vector field T a is orthogonal to 
the foliation 

The strategy of the proof is the same as usual: the geometric objects are split over the basis 
(/°> e t)i=i,2.3 where f a is transverse to S t and (e°)i=i.2,3 tangent to S t . 
Considering the 4-form 

/„ A e? A e a 2 A e a 3 

f a is decomposed as follows: 



r = (3e%- £ 5'et 



i=l,2,3 



with 

S' 1 = 9cdf c ef and /3 2 = 1 + ]T (S 1 ) 2 



i=l,2,3 



The volume form fi[g] satisfies: 



M = -oh A e? A el A e% 



and its contraction with the vector f a is: 



f a ^[g] =\(3e a - I Mtt = Pel Ae a 2 Ae a 3 

i=l,2, 3 
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as a consequence, we have: 

i%(*f a ) = \\f*\\%%(*f a ) = \\f a \\f a M$ = l|T a ||/3 e ? A e a 2 A eg. 
We then deal with V c </> which can be written: 

V c 4> = bf - o'ef 

1=1,2,3 

where 

, 9abV^4 / Q Va^ + E. t= l,2.3^<Va^ , .A, „, . 

o = = and a\ = e^\/ a <p — bo tor i g {1, 2, 3}. 

P P 

Consequently, its norm is: 

V 0V C = 2^ ( e i V a0) 

" i=l,2,3 



and the restriction to S# of -kVh<f> is: 



/? 2 



.* / tt i\ w ^1=1,2,3 "i ' .* / j. s 

4 t (*V b 0) = 4 t (*/a) 



i=l,2,3 

Using these results, the energy 3-form is given by 



rv a 0+ J2 ^fV^| e?A | Aeg . (A.7) 



1 / (/"Va^+Ei^XS^efVa 



e? A e£ A eg 



i=l,2,3 



E ( e i^) 2 +T +6 T ) H!T-I|/3ef Aeg Aeg 



2 I ^2 ^ v-i • »r/ I ■ 2 ■ - 4 

r i=l,2,3 



i=l,2,3 

To get the equivalence, it is sufficient to notice that: 



2 A&\ \\rra\\1 



2 4 / /3 



-e? A eg A eg 



• gcd,T c f d is a positive function over the compact V and, as such, is bounded below and 
behind by two positive constants; 

• as already noticed in remark f2.1 91 j3 is a positive function over V; 

• and, finally, the scalar products j3, <5i and the difference f3 2 — Si are bounded below by 1 
and above by a certain constant since we are working on a compact set. 

The energy equivalence then holds: 

/ <h, M - / + i £ «v„tf + f. + ^ 



2 4 

i=1.2,3 



□ 
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i+ 




{t = 0} 



S = {t = e} 



A. 3 Local characteristic Cauchy problem for small data 

It has been considered as interesting to give here a proof of the existence of a solution for small 
data because of its technical simplicity. 



neither i° nor i + . Let r be a "reverse" time function on M, in the sense that its gradient is past 
directed with respect to T a (defined in section l2.2.2p . We assume that r(i + ) =0. M is endowed 
with an orthonormal basis (e?)j = o. 1,2.3 such that eg is colinear to Vt and (e?)^!^^ is tangent 
to the time slices {t = constant}. The integral flow associated with d T is denoted by $ T . 

Let t be a positive constant smaller than 1 and consider <f> e (z + ). Let S be a uniformly 
spacelike hypersurface for the metric g between J? + and {$ £ (p)|p £ We assume that § is 

uniformly spacelike, transverse to ^ + in the past of the support of the characteristic data, 9, 
and contains $ £ (* + )- 

Remark A. 8. The geometric framework is then exactly the same as in section Wlk the timelike 
vector field is not normal to the foliation. Nonetheless, since the hypersurface S is uniformly 
spacelike, the same estimates as in section [2.3\ hold without the nonlinearity. 

Finally, the future of S in M is foliated by the surfaces § T = {Q e - T (p)\p <G S} for r in [0, e] 
so that § = S e and §0 = {i + }- The future of § is denoted by "R and the subset of M between §0 
and § r , D? r . 

The solution of the nonlinear problem is approximated via solutions of the linear problem 
on . Hormander solved this problem in |32) : 

Proposition A. 9 (Hormander). Let us consider the linear inhomogeneous characteristic Cauchy 
problem on ,y + : 



where f is in Lj^H 1 and 9 is a function whose compact support does not contain i + or i . Then, 
this problem admits a unique global solution in the future of £0 in C°([0, e], H 1 (§ T )). 

Using proposition IA.9I and estimates for the linear problem, the following theorem holds: 

Theorem A. 10. Let us consider the nonlinear characteristic Cauchy problem on J^ + ; 





□0 + 
^1 



gScal§(f> = f 
= 9e H\J+). 




Ou + ^ScalgU + bu 3 = 
uL + =9^H\J+). 
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where is a function whose compact support does not contain i + or i . Then, for \\0\\h 1 (J'+) 
small enough, there exists a uniformly spacelike hypersurface § close enough to + such that 
this problem admits a smooth global solution on 31 in C°([0, e], H 1 (S T )). 

Remark A. 11. The proof of the well-posedness in C°([0, e], H 1 (§ T )) is given in section \3.1.S\ 
were the geometric estimates required to obtain it are established (see theorem \3.5\ which remains 
true in that context). 

Proof. Let uq be a solution on 31 of the problem: 

U4> + iScalg</> = 

and let (u n ) n eN be the sequence of smooth functions on 31 defined by the recursion: 

□ti„+l + iScalgU„ + i + bul = 

««+iU =6eH l {y+). 

The sequence defined by the difference of two consecutive terms of this sequence is denoted by 
(S n = u n+ i — u„)„ s n- For neN, the smooth function S n satisfies the Cauchy problem: 

OS n + iScalg(5„ = -b(ul + u„it„_i + i4-i)<5»-i 
S n = on S+. 

The proof of the convergence is made in two steps: the first one consists in proving that, for 
initial data which are small enough, the sequence (u n ) is bounded; the second part proves the 
convergence of (u n ) by showing that the sequence (S n ) is summable. 

Proposition A. 12. For \\0\\H 1 (y+) small enough, the sequence swp T£ ^ e i ||it n ||/fi(s T ) is bounded. 

Proof. Let n be a integer greater than 1. Let us finally consider the energy tensor associated 
with the linear wave equation: 

Tab = VaWn+lV&Mn+l + g a b f~ ^ V c "n+1 V C U„+i + -y^- 

The energy associated with a time slice S T is written as: 



EuM= f I J2 ^aU n f + ^ I d/i, 
T \i=0,l,2,3 J 



and it is equivalent to J„ i*(*egT a ;,) (with constants which only depend on the geometric data, 
the L°°-norm of b and the Killing form of eg). 

The error term associated to this energy tensor is: 

V a (e b T ab ) = W {a e b Q ] T ab + u n+1 e^V a u n+1 - -Scal$ti n+ iegV Un+i - (egV a u n+1 )bul, 
which is smaller than, in absolute value: 

\V a e b T ab \ < C ^>?V a n n +i) 2 + u 2 n+1 ] + Cbu 6 n , 

where C is a positive constant depending on sup (|Scal§|), sup M|V^ a eQ^||^ and the foliation S r . 

The next step consists in using the Sobolev embedding of H 1 (§ T ) in L 6 (S T ). There exist two 
obstacles to the use of this embedding: 
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a. the first is the fact that the estimates of the nonlinearity must be uniform over the foliation 
in the sense that it must not depend on the parameter r of the foliation (or the Sobolev 
constant must be the same all along the foliation) and does not explode despite the behavior 
of the Sobolev constant ; 

b. the second comes from the fact that we must deal with the singularity in i + . 

To deal with the second problem, the manifold M is extended beyond J 2 " 4 " by pulling back- 
wards the hypersurface S through the flow associated with the vector field d T of the time function 
t. Since the regularity of the metric is arbitrarily smooth at i + (say, at least C 2 , in order to 
insure the existence of the different curvatures), this gives an extension as a smooth Lorentzian 
manifold of the manifold (M, g) in the neighborhood of i + . J> + is then the past light cone from 
i + obtained from a C 2 -extension of the metric g behind i + . 

To obtain a Sobolev embedding from H 1 (§ T ) into L 6 (S T ) uniformly in r, it is necessary to 
control the blow-up of the Sobolev norm. This is achieved in the same way as in section I3.i.2l 
using proposition ^. 31 A direct consequence of the decay of the function b as stated in assumption 
|d]is then the following inequality: 

\W a e b T ab \^ T < CE Un+1 (§ T )+C\\u n \\ 6 HH§T) 

< CE Un+1 {§ T ) + C I sup E Un (S T )) . 

\re[0,e] J 

As a consequence, there exists a constant C which depends only on the scalar curvature, the 
Killing form of eg, the decay of b and Sobolev constants such that: 

E Un+1 (Sr) <C\ f E Un+1 (§ r )dr + \\9\\ 2 H1{ ^ +) + I sup (E Un (S r )) 
wo yreltv 1 

Remark A. 13. The constant C can be chosen arbitrarily high. As a consequence, it is rescaled 
later without consequence for the proof (see remark \A.15\ in the proof of proposition \A.lJ$ . 

Using Stokes theorem and Gronwall lemma, as in section \2. 2.21 the energy of satisfies: 
E Un+1 (S T ) < Cexp(Ce) ^\\0\\ 2 m ^ +) + sup^ u „(S T ) X 
For n = 0, we have: 

E U0 (S T )<CcMCe)\\e\\ 2 Hl ^ +) . 
We denote by (C„) the sequence defined by: 

C n = sup {E Un (S T )} . 

r£[0,c] 

This sequence satisfies the inequality: 

Vn 6 N,C n+1 < Cexp(C'e) C 3 n + \\6\\ 2 Hl{y+) | with C < Ccxp(Ce) \\e\\ 2 Hl{ ^ +) . 

° \ ~y 

Let us then consider the sequence (c n )„ defined by: 

Co = afJ 
Cn+i = a (cl + (3) 



52 



The purpose is to choose correctly such that this sequence is bounded. The function 



x i ^ a [x 



~X + satisfies: 



/3) has three fixed points provided that the discriminant of the polynomial X 3 



/3 2 



27a 3 ~ 27 
Using Cardano's formulae, its three zeros are: 



<0, ie->((7exp(C6)) 3 ||0|| 4 ffl( ^ +) . 



An 



Ai = 



f 



3a 
3a 

T 



cos — arccos 



cos 



cos — arccos 



arccos — 



27/3 2 a 3 



27/3 2 a 3 



27/3 2 a 3 



+ 



2tt 

47T 

T 



(A.8) 

(A.9) 
(A.10) 
(A.ll) 



Since 



27/3 2 a 3 



- arccos 



arccos 



27,3 2 a 3 



27/3 2 a 3 



2- 



— € 



-in 



7T 7T 

6' 3 
6 ' 



37T 57T 



these roots can be compared as follows: 

Ai < < A 2 < A . 

The fixed points An and Ai are repulsive whereas the fixed point A2 is attractive. As a con- 
sequence, if the (positive) initial condition en is below the positive repulsive fixed point (the 
greater fixed point An) of the function x *- > a (x 3 + /3), that is to say if 



cos — arccos 



27^ 2 a 3 



3 cos 



1 / /27^ 2 a 3 

3 arccos I -y — 4 — 



> aj3 



27/3 2 a 3 



the sequence (c„) converges towards A2. This inequality is always satisfied as soon as: 

4 



27a 3 



< 0. 



As a consequence, assuming that 



27(C*exp(Ce)) 3 



> 



the sequence (c„) converges to the remaining attractive fixed point A2; (c„) is bounded and so 
is (C„), which is the expected result. 
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Another useful consequence of the convergence of the sequence (c„) is the following. The 
limit of this sequence satisfies: 



A 2 < 



4 



3a y 3Cexp(Ce) 
As a consequence, there exists a integer uq such that: 



Vn > n , sup {E U JS T )} < c n < 2j 4 (A.12) 
rG[o,e] y 36 exp(Ge) 

□ 

Proposition A. 14. The sequence (u n ) converges on 3? in C°([0, e], H 1 (§ T )), that is to say in 
the norm ^sup re [ c j |M|^-i( S j 

Proof. The method is exactly the same as in the previous proposition. Let n be a positive integer 
and consider the energy tensor associated with the linear wave equation for the function S n 

( 1 51 
Tab = V a <5 n V 6 <5„ + g a b -t:V c (5„V c (5, 



2 " 2 

The energy associated with a time slice § T is written as in the previous proposition: 



EsM= [ \\ J2 (e?V<5„) 2 + f ]d^ 

T \ »=0, 1,2,3 J 



and it is equivalent to Jg «g (*egT a t,) (with constant which only depends on the geometric data, 
b and the Killing form of eg). 
Finally, the error term is: 

V a (e b T ab ) = V {a e$T ah + S n e%X7 a 6 n - ^Scal~ g S n e%V a 6 n - &(egV„<5 n )<5 n _i(^ + u„u„_i + u^). 
and can be estimated in absolute value by: 

\V a e b T ab \^ T < CE Sn (S T ) + 2 f + K-i)^- 

where C is a positive constant depending on sup (|Scal§|), sup ^| |V^ a e^ | \j , sup(|5|) and the 
foliation S r . 

Using Holder inequality, proposition 13.31 for the foliation § T , the non-linearity in the error 
term is estimated by: 

2 12 

b2 $l-i u tiVSr < (supb) (/ 5j l _ifi§ T J (/ u 6 n fi Sl 
s T \s T J \Js T / \Js T 



< c 



(sup Sr b) 2 



T 6 



ll<Mlff 1 (S T )ll' lt ™|lff 1 (S T ) 



< C||<5Ti|| fl -i(s T )ll u n||_H-i(S T ) 

The same inequality holds for 6^ l _ 1 uf l _ 1 . 
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Finally, there exists a constant K, such that: 

/ \V a e b Q T ab \fi ST 

<K\ [ E Sn (§r)dr + el sup (Si B _,(S r ))] sup ( sup (E Un (S r ))\ 

WO \re[0,e] / k>n-l \r£[0,e] J 

Stokes theorem is then applied beween § r and J^ + : since the characteristic data for S n are zero, 
the only remaining term is the energy on the surface § r . Modulo a constant which only depends 
on the same data as the constant C, the integral inequality holds: 

E Sn (S T )<K\ [ E Sn (B r )dr + el sup (£*„_,(§,-))) sup ( sup (E Ufl (S r )) 

WO \rG[0,e] J k>n-l \re[0,e] y 

for some contant K and, using Gronwall's lemma, we get: 

EsM < K exp(^e)e ( sup [E &n _, (S r )) ) sup f sup (E Un (S r ))) . 

\re[0,e] / k>n-l \re[0,e] J 

Remark A. 15. The constant K, as the constant C depends only the foliation S T , its scalar 
curvature, the Killing form of eg and the supremum of b. As a consequence, up to a rescaling 
of C or K , these constants can be chosen to be equal. 

Finally, the sequence (S n ) satisfies: 

sup (E Sn (S T )) < C exp((7e)e ( sup (S 5n _ 1 (S r )) j sup ( sup (E Un (S r ))) . 

re[0,e] \r-e[0,e] J fc>n-l \re[0,e] J 

Using inequality (|A.12[) , we have: 



Cexp(Ce)e sup sup (E Un (S r )) ) < Ceexp(Ce) 2j-~ 

k>n-i \re[o,e] J \ y3Gexp(Ge, 

< ye. (A.13) 

Since e is smaller than j^, the sequence (S n ) is then eventually contracting. The series of 
(sup,, e [ e ] (||<5„||ffi(s r )) 2 )n converges in the norm ^sup T6 [ e j IMIffi( §r )) , that is to say in 
C°([0, e], _ff 1 (S T )), and so does the sequence (u n ). □ 

End of the proof of theorem ] A. 1 (A The proof of the local existence is a direct consequence of 
the fact the sequence (u n ) n ^jq converges strongly on 3? for the norm 

SU P (IMIffMs,-)) ■ 

re[0,e] 

Let u be the limit of the sequence (u n ) n . The only remaining thing to show is that the limit 
solves the problem of theorem lA.101 It is clear that u satisfies the initial conditions since all the 
functions u n are identically equal to 9 on J? + . Finally, when noticing that u is in 7J 1 ([R) which 
is continuously embedded in L 3 (D?) (since H is four dimensional), the sequence (u^) n converges 
in L 1 ^) and, as a consequence, in the distribution sense, u then satisfied the equation 

1 , 

Du + - Seals u + bu A = 
6 

in the distribution sense. □ 
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